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Preface 


The  purpose  of  the  thesis  was  to  qualitatively  describe  the  effects  of  aberra¬ 
tions  on  off-axis  beam  steering  for  multi-aperture  systems.  A  new  expression 
analogous  to  Goodman’s  thin  lens  was  developed  which  shows  the  beam  degrada¬ 
tion  that  occurs  at  angles,  with  respect  to  the  optic  axis,  much  greater  than  five 
degrees  for  a  thin  lens,  and  two  degrees  for  a  thick  lens.  It  is  my  hope  that  this 
thesis  will  be  the  first  in  a  series  which  will  not  only  prove  and  validate  the  idea 
of  phased  optical  beam  steering,  but  to  make  an  engineering  analysis  on  such  a 
system. 

The  analysis  described  in  this  research  is  of  special  interest  to  organizations 
such  as  SDI  since  it  defines  regions  and  specifications  for  a  multi-aperture  system 
which  will  not  work.  Hopefully  with  this  research,  a  reliable  system  can  be  made 
to  achieve  the  goals  of  SDI. 

I  would  like  to  thank  my  advisor,  LtCol  Jim  Mills,  who  has  spent  many  an 
hour  discussing  this  project  with  me.  I  would  also  like  to  especially  thank  Maj. 
Jim  Lupo  (user  friendly)  for  his  invaluable  support,  as  confidante,  and  computer 
expert  -  without  whom  I  am  sure  there  ..ould  have  been  a  baseball  bat  through 
more  than  one  terminal.  I  would  like  to  thank  my  parents,  who  happened  to  be 
stationed  here  three  days  after  I  did  -  please  don’t  follow  me  to  my  next  assign¬ 
ment!  Its  awful  hard  to  have  a  daughter  around  for  the  first  time  in  five  years, 
and  to  refrain  from  calling  because  of  thesis.  And  most  importantly,  I  can  never 
thank  these  two  people  so  that  they  will  ever  understand.  To  the  person  who 


stood  by  me  to  help  me  finish  this  thesis  -  even  to  staying  up  all  night  with  me 
the  night  before  it  was  due.  May  our  future  together  be  as  beautiful  as  our  past. 
And  to  my  bestest  friend  and  teddy  bear,  I  can  only  hope  that  I  helped  you  as 
much  as  you  helped  me  as  we  went  through  this. 


Christina  N.  Walton 
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Abstract 


The  effect  of  aberrations  on  off-axis  beam  steering  for  multi-aperture  systems 
was  qualitatively  investigated.  A  new  expression  analogous  to  Goodman’s  Thin 
Lens  is  proposed  which  shows  the  beam  degradation  that  occurs  at  angles,  with 
respect  to  the  optic  axis,  much  greater  than  five  degrees  for  a  thin  lens,  and  two 
degrees  for  a  thick  lens.  This  expression  was  then  applied  to  three  multi-aperture 
configurations  as  a  way  to  correctly  predict  the  field-of-view  for  various  optical 
systems.  Analytical  predictions  supported  by  computational  and  experimental 
results  are  included.  Although  the  beam  degradation  as  the  field-of-view  or  field 
angle  is  increased  has  a  dramatic  effect  on  the  image  quality,  the  power  loss  as  the 
field  angle  is  increased  will  be  the  major  factor  in  designing  a  system  of  this  type. 
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Aberrations:  Their  Effect  On  Off- Axis  Beam  Steering 


I.  Introduction 


1.1.  Problem  Statement 

In  this  thesis,  the  effect  of  aberrations  on  off-axis  beam  steering  for  multi¬ 
aperture  systems  was  qualitatively  investigated.  A  new  expression  analogous  to 
Goodman’s  thin  lens  (1:80)  is  proposed  which  shows  the  beam  degradation  that 
occurs  at  angles,  with  respect  to  the  optic  axis,  much  greater  than  five  degrees  for 
a  thin  lens,  and  two  degrees  for  a  thick  lens.  This  expression  was  then  applied  to 
three  multi-aperture  configurations  as  a  way  to  correctly  predict  the  field-of-view 
for  various  optical  systems.  This  thesis  includes  analytical  predictions  supported 
by  computational  and  experimental  results. 


1.2.  Air  Force  Applications 

Today  a  single  phased  array  radar  can  do  what  previously  might  have 
required  a  battery  of  mechanically  steered  dishes.  Phased  array  radars  are  a  series 
of  point  emitters  which  radiate  at  microwave  frequencies.  If  a  phase  shift  is  intro¬ 
duced  at  each  emitter,  then  the  beam  can  be  steered  as  opposed  to  rotating  a  plat¬ 
form  (or  radar  dish).  The  replacement  of  a  single,  movable  antenna  with  an 
immobile  array  of  radiating  elements  can  offer  more  advantages  than  just  elec¬ 
tronic  steering,  reliability  being  one  of  them.  A  fixed  array  does  not  depend  on 
fallible  mechanical  components  such  as  bearings  and  motors.  With  increasing 


simplicity  and  compactness,  concomitant  reductions  in  weight  and  increases  in 
reliability,  phased  array  radars  can  be  deployed  in  space  within  a  few  years  (2:9-4- 
102). 

The  President's  Strategic  Defense  Initiative  (SDI)  has  been  proposed  as  the 
central  defense  concept  in  our  nation  for  the  years  to  come.  One  of  the  SDI  mis¬ 
sions.  midcourse  discrimination,  requires  simultaneous  acquisition,  identification, 
and  tracking  of  a  large  number  of  targets  over  a  twenty  degree  field-of-view.  The 
obvious  benefits  of  a  phased  array  multi-aperture  optical  system,  therefore, 
become  very  important  to  the  Air  Force  and  the  Department  of  Defense. 

The  ability  to  achieve  a  large  field-of-view  from  a  multi-aperture  optical  sys¬ 
tem  has  been  doubted  by  many  based  on  the  transfer  function  for  a  thin  lens 
derived  by  Goodman  (1:77-80).  Based  on  the  new  expression  derived  to  include 
the  effects  of  off-axis  radiation,  beam  steering  for  an  uncorrected  multi-aperture 
optical  system  is  not  realizable. 


1.3.  Background 

The  transmission  function  for  a  thin  lens  is  given  by  Goodman  (1:80)  as: 


f,(jr.y)=exp 


-jk 


1  o  o.. 

D'+r) 


(i.i) 


The  expression  was  derived  to  enable  diffraction  analysis  for  on-axis  radia¬ 
tion.  This  function  was  used  to  show  how  a  lens  can  perform  the  Fourier 


transform  of  an  input. 


Geary  and  Peterson  (3:52-54)  derived  a  new  transfer  function  for  a  thin  lens 
which  included  one  more  term,  from  the  binomial  expansion,  that  had  been 
assumed  negligible  in  Goodman's  expression.  They  showed  that  this  extra  term 
included  the  effects  of  spherical  aberration.  However,  if  the  rest  of  the  terms  are 
kept  in  the  expression  for  the  transfer  function,  no  further  information  is  gained 
about  the  other  third  order  aberrations  such  as  astigmatism,  coma.  Petzval  field 
curvature,  and  distortion.  The  only  information  found  is  an  exact  specification  of 
spherical  aberration. 

In  none  of  the  material  found  to  date  is  an  exact  specification  for  a  transmis¬ 
sion  function  of  a  lens  which  is  appropriate  for  off-axis  radiation  (further  than  the 
paraxial  approximation).  The  transmission  function  for  the  lens  is  a  function  of 
the  distance  the  light  travels  through  the  lens,  relating  that  distance  to  a  phase 
lag  the  light  undergoes.  Light  incident  on  a  lens  surface  at  an  angle  larger  than 
zero  degrees  (off-axis),  will  travel  a  further  distance  than  was  previously  assumed. 
This  additional  distance  actually  yields  far  larger  aberrations  than  were  previously 
expected  for  the  optical  system,  thereby  significantly  degrading  system  perfor¬ 
mance. 

1.4.  Multiple  Apertures 

A  multiple  (or  synthetic)  aperture  is  formed  when  separate  optical  systems 
are  combined  to  function  as  a  single  larger  aperture  (4:2).  Coherent  imaging 
through  a  synthetic  aperture  has  many  advantages,  among  them  ease  of  manufac¬ 
turing  and  reliability  (4:2-4).  The  more  apertures  contained  within  a  multi- 
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aperture  system  the  greater  the  light  gathering  capability.  The  central  peak  of  a 
multi-aperture  system  is  more  narrow  than  that  of  a  single  optic  which  has  the 
same  outer,  or  effective  radius.  This  is  advantageous  to  tracking  or  two  point 
resolution  as  was  pointed  out  by  Watson  (5:40-62). 

Coherent  imaging  through  a  synthetic  aperture,  where  the  dc  component  of 
the  light  has  been  blocked,  also  has  the  added  feature  of  being  an  edge  enhancer. 
Therefore,  there  will  be  "ringing"  around  any  edge  of  an  object  coherently  imaged 
through  an  optical  system  as  shown  in  Figure  1.  Bergey  (6:72-75)  pointed  out 
that  these  patterns  could  be  controlled  by  the  number  of  apertures  in  the  optical 
system  as  well  as  the  placement  of  the  apertures  in  the  system. 


aberrations:  spherical,  astigmatism,  coma,  distortion,  and  Petzval  field  curvature, 
will  be  discussed  in  this  thesis.  Spherical  aberration  is  dependent  upon  the  fact 
that  the  lens  is  spherical  in  shape  as  opposed  to  aspheric  (e.g.  hyperbolic).  .Astig¬ 
matism  and  coma  are  off-axis  effects.  If  the  incident  angle  of  the  plane  wave  with 
the  optic  axis  is  zero,  then  the  effect  due  to  these  two  aberrations  is  zero.  .As  the 
incident  angle  grows,  the  effects  due  to  these  aberrations  grow  as  well.  This 
greatly  affects  the  final  image  since  the  greater  the  aberration,  the  greater  the 
beam  deformation.  Such  deformations  cause  a  decrease  in  field-of-view.  Field 
curvature  and  distortion  only  affect  the  image  of  an  extended  object.  Since  this 
thesis  deals  only  with  the  impulse  response  of  the  optical  systems,  these  two  com¬ 
ponents  of  the  third  order  aberrations  will  be  neglected. 

When  the  affects  of  these  aberrations  are  included  in  a  multi-aperture  system, 
the  beam  degradation  as  the  incident  angle  grows  will  become  even  more  of  a 
problem.  The  problem  with  two  point  resolution  in  multi-aperture  systems  is  the 
addition  of  side  lobes,  which  can  be  as  large  if  not  larger  than  the  central  peak 
(5:40-62).  .As  the  incident  angle  grows,  the  lobes  on  one  side  of  the  intensity  pat¬ 
tern  grow  much  faster  than  the  lobes  on  the  opposite  side.  This  increases  the  pos- 
siblity  of  error  in  reading  the  pattern.  Therefore,  the  beam  deformation  problem 
becomes  worse  in  a  multi-aperture  system.  Aberrations  are  one  major  limitation 
to  the  field-of-view  for  a  multi-aperture  optical  system. 
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1.6.  Scope 


For  SDI  purposes,  optical  phased  array  radars  are  an  optimal  solution  lor 
tracking  and  targeting.  The  requirement  of  a  twenty  degree  field-of-view  means 
that  for  far-field  imaging  a  plane  wave  could  be  impinging  the  lens  at  an  angle 
with  respect  to  the  optic  axis  of  twenty  degrees  or  more.  The  transmittance  func¬ 
tion  for  the  Goodman  thin  lens,  which  has  been  assumed  in  all  previous  work  is 
not  applicable  for  angles  in  excess  of  two  degrees,  as  the  experimental  results 
showed.  Goodman,  in  defining  his  thin  lens  assumed  the  paraxial  approximation. 
In  this  thesis,  a  new  definition  for  the  thin  lens  will  be  derived  which  includes  the 
angle  of  incidence  implicitly.  The  only  approximation  in  this  method  was  to 
assume  that  the  angle  change  at  the  surface  due  to  the  index  of  refraction  is  negli¬ 
gible. 

This  thin  lens  expression  was  applied  to  three  multi-aperture  systems  to  show 
how  the  beam  degradation  affects  the  field-of-view  as  the  incident  angle  is 
increased.  The  three  multi-aperture  systems  used  are  shown  in  Figure  2.  Two  of 
the  three  multi-aperture  systems  are  being  tested  at  the  Air  Force  Weapons 
Laboratory,  Ivirtland  Air  Force  Base,  New  Mexico.  The  three  aperture  system  is 
called  Phasar  (7)  and  has  been  in  operation  for  approximately  two  years,  while 
the  four  aperture  system  is  called  the  Multiple  Mirror  Tracking  Telescope 
(MMTT)  (7).  which  will  begin  operation  in  January  1988.  The  six  aperture  sys¬ 
tem  is  in  the  same  configuration  as  the  Multiple  Mirror  Telescope  (MMT)  in 
Arizona,  where  it  is  engaged  in  astronomical  and  other  scientific  endeavors.  Bor- 
gey  showed  that  the  six  aperture  system  is  one  of  the  optimum  designs  for  multi- 
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aperture  systems  (6:72-75).  Also,  Watson  (5:40-62)  demonstrated  that  six  aper¬ 
tures  has  a  better  two  point  resolution  than  the  other  systems  he  investigated. 
All  of  the  systems  used  in  this  thesis  are  symmetric  with  the  optical  axis  as  shown 
and  the  lenses  are  separated  by  a  small  distance  to  avoid  overlapping  in  the  com¬ 
putational  analysis  as  well  as  to  allow  for  the  fixed  supports  these  systems  would 
require  in  reality. 
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This  thesis  will  deal  strictly  with  coherent  imaging.  In  Chapter  2.  aberra¬ 
tions  and  Goodman’s  thin  lens  will  be  discussed.  The  new  lens  transmittance 
function  will  be  derived  and  applied  to  two  different  lens  designs  in  Chapter  3. 
Chapter  4  will  verify  the  methods  applied  by  comparing  the  previous  results  with 
the  geometrical  rav  trace  program  FALCON  and  in  Chapter  5.  experimental 
analysis  is  used  to  support  the  work  done.  Chapter  6  will  apply  thesis  lens 
designs  to  the  three  multi-aperture  systems.  A  comparison  of  these  systems  will 
be  made  in  Chapter  7  and  the  conclusions  drawn  will  be  presented  in  Chapter  S. 


II.  Aberrations  and  the  Perfect  Thin  Lens 


This  chapter  develops  the  concept  of  aberrations  and  the  perfect  thin  lens. 
Aberrations,  as  discussed  by  Born  and  Wolf  (8:459-490),  will  be  applied  to  an  opt¬ 
ical  model.  The  aberration  coefficients  are  then  derived  for  practical  comparisons 
of  these  optical  models,  and  the  case  of  Goodman’s  thin  lens  applied.  This  model 
was  computationally  analyzed  and  the  resulting  intensity  patterns  are  given  at  the 
end  of  this  chapter.  In  addition,  the  lens  parameters  of  the  thin  and  thick  lenses 
used  in  all  of  the  analyses  in  this  thesis  are  included. 


II. 1.  Aberrations 

The  following  development  follows  closely  that  of  Born  and  Wolf  (8:460-472). 
Consider  the  centered  optical  system  shown  in  Figure  3: 

In  order  to  get  a  perfect  Airy  pattern  (8:460)  in  the  image  plane,  a  perfect 
but  truncated  spherical  wave  must  exit  the  lens  system  (at  the  exit  pupil).  How¬ 
ever,  if  a  wavefront  exits  the  lens  system  and  is  different  from  a  perfect  spherical 
wave,  then  aberrations  exist  which  are  of  the  form  Q  -  Q,  (the  distance  between 
the  two  wavefronts).  It  is  convenient  to  regard  6  as  the  total  effect  due  to  aber¬ 
rations.  where  6  is  of  the  form: 

d  =  6(Y;.P.e)  (2.1) 

and: 


V'  =  object  height  from  the  optic  axis  in  the  image  plane 
p  —  radial  distance  from  the  optic  axis  in  the  plane  of  the  exit  pupil 
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Gaussian 
Reference  Sphere 


Wavefront 


Figure  3.  Reference  Frame  and  Notation  for 
Aberration  Analysis 

9  =  angle  of  rotation  around  the  optic  axis  in  the  plane  of 

the  exit  pupil 

Including  aberrations,  the  Huygens-Fresnel  equation  becomes: 

OU„„)-  UL(-,fc)//(HEa±tlils  (2.2) 

A  Z  S 

Where  dS  is  the  spatial  integration  over  the  exit  pupil,  lambda  is  the  wavelength 
of  interest,  k  is  the  wave  number,  and  z  and  s  are  optical  path  differences  from 
the  object  to  image  planes. 

The  function  d  can  be  separated  in  many  fashions.  Two  of  the  most 
famous  ways  are  given  by  Seidel  and  Zernike.  The  Seidel  aberrations  are  defined 
for  geometrical  optics,  whereas,  the  Zernike  circle  polynomials  were  defined  for 


diffraction  analysis.  The  representations  for  the  primary  aberrations  in  both 
forms  is  given  below  in  Table  II- 1,  where  the  expansions  for  6  are  given  by: 


Seidel-.  <t>  =  At’nm  pn  cosm  9  (2.3) 

Zernike:  6  =  enmAlnmR?(p)cosrn9  (2.4) 


TABLE  n-1. 

Representation  of  the  Primary  Aberrations 


Aberration 

1 

n 

m 

Seidel  Form 

Zemike  Form 

Spherical 

0 

D 

0 

^040  P4 

-j=-AmR04(p)  =  ~AM0(6p^  -  6p2  +  1) 

Coma 

0 

3 

1 

Ao3ip3cos0 

A  031#3  (p)COS0  =  A  031  (3p2  -  2p)cos0 

Astigmatism 

0 

2 

2 

Aq22  p2COS29 

AO22^2(p)cos20  =  A022p2(2cos29  -  1) 

Field  Curvature 

1 

2 

0 

^120P2 

•^l»*°(P>  =  ^U0<2pJ-l) 

Distortion 

1 

1 

1 

A  in  pcos0 

A 1 1 1 R  i  (p)cos0  =  A  i !  i  pcosO 

From  Eq.  (2.2),  the  transfer  function  G  can  be  defined  as  follows: 

G(*o-yo)  =  exp  (iko)  (2.b) 

where  G  contains  all  of  the  contributions  to  the  image  due  to  aberrations.  There¬ 
fore,  all  of  the  aberration  effects  produced  by  an  optical  system  can  be  combined 


said  that  the  optical  system  is  perfect,  but  there  is  a  transparency  across  the  exit 
pupil  .  G,  which  contains  all  of  the  aberrations.  The  sole  effect  of  these  aberra¬ 
tions  is  seen  to  be  the  introduction  of  phase  distortions  within  the  bandpass  of  the 
system.  These  phase  distortions  can  have  a  severe  effect  on  the  fidelity  of  the  sys¬ 
tem  (1:121). 


II.2.  Lens  Design 

Two  lenses  were  used  in  the  experimental  analysis.  These  lenses  were  then 
used  to  test  out  the  analytical  and  computational  portions  of  this  thesis.  The 
first  lens  was  a  thin  lens  which  was  chosen  to  verify  the  computational  analyses, 
while  the  second  lens  was  chosen  as  a  thick  lens  to  show  the  limitations  of  the 
analyses.  Given: 

LENS  1: 

/?!  =  16.84cm 
R  2  =  — 44.23cm 
H  =  4.372cm 

On— axis  thickness  =  ,794cm 

and 

n  =  1.5 

Where  n  is  the  index  of  refraction  for  any  non-specific  glass.  The  second  lens  was 
designed  to  test  the  limits  of  the  approximations  made.  Lens  2  is  a  thick  lens 


which  will  have  large  aberrations  present: 


I 


LENS  2: 


and 


R  i  =  11.37  cm 
R2  =  —  11.63cm 
H  =  6.0976cm 
On  —  axis  thickness  =  3.5 cm 


n  =  1.5 

The  aperture  stop  in  both  optical  systems  is  1  cm  in  radius. 


II. 3.  Perfect  Thin  Lens 

A  lens  is  said  to  be  a  thin  lens  if  a  ray  entering  at  coordinates  (x,y)  on  one 
face  emerges  at  approximately  the  same  coordinates  on  the  opposite  face,  i.e.,  if 
there  is  negligible  translation  of  the  ray  within  the  lens  (1:77).  Figure  4  gives 
the  description  of  the  lens  used  in  this  discussion.  This  discussion  follows  Good¬ 
man  (1:77-85). 
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Figure  4.  Goodman's  Thin  Lens 

The  phase  delay  given  to  an  incident  wave  going  through  this  lens  is: 

6(x,y)  =  kn^x,y)+k(V0-V(x.y))  (2.6) 


Where  n  is  the  index  of  refraction,  the  first  term  is  the  phase  delay  introduced  by 
the  lens  and  the  second  term  is  the  phase  delay  introduced  by  the  remaining 
region  of  free  space  between  the  two  planes.  Therefore,  the  transmittance  func¬ 
tion  of  the  thin  lens  is  given  by: 

t,(x,y)  =  exp(j*V0)exp(;7;(n— l)V(x,j/))  (2.7) 

And,  the  complex  field  immediately  behind  the  lens  is: 

Ut'{x,y)  =  t,{x,y)U,{x,y)  (2.8) 

Thus,  the  problem  is  to  define  the  thickness  function  of  the  lens,  'T(x.y)  .  To 

define  this,  break  up  the  lens  into  two  separate  parts,  as  shown  in  Figure  5.  Then 

the  problem  is  to  solve  for  each  half  of  the  distance,  or  phase  delay,  that  the  lens 

imparts. 

After  solving  for  distances  the  light  travels  in  the  lens,  and  V2  are  given 
by: 

V,(*,y)  =  v01-(i? ,-(/?  r  -X2~y2]  (2.9) 

_1_ 

V2(x,y)  =  V02-(-/?o-[/?| -x2-y2]  2 )  (2.10) 

Adding  the  terms  together,  remembering  that  V01+V02=V0: 


V(x,y)  =  V0-/?i(l—  [i_i£^!l^)+/?2(i_  (2.11) 

The  most  important  approximation  taken  in  Goodman's  analysis  is  the 
paraxial  approximation.  After  completing  the  binomial  expansion  of  the  two 
square  roots,  this  approximation  requires  that  only  the  first  two  terms  are 
retained  and  all  other  terms  are  ignored.  This  assumes  that  the  radii  of  the  two 


spheres  which  make  up  the  thin  lens  are  much  greater  than  the  radius  of  the  lens 


itself  (or  the  distance  from  the  optic  axis).  After  this  assumption  is  made,  the 


Figure  5.  Calculation  of  the  Thickness  Function 
new  equation  for  V  is: 


Hr  ,y) 


(Hfy2) 

f  i 

1 

o 

F, 

r2 

(2.12) 


Substituting  this  expression  into  Eq.  (2.7)  and  remembering  the  len's  makers 


equation: 


The  final  equation  for  the  perfect  thin  lens  is  then: 


=  exp  ~-(x"+r) 


(-•1-1) 


Where  the  exponential  phase  factor  eip(jknV 0)  can  be  ignored  since  we  are  only 


looking  at  the  intensity  patterns  in  this  thesis.  When  this  equation  is  used  to  find 


the  complex  field  amplitude,  (ignoring  the  phase  factors  in  front  of  the  integral 


since  they  will  negate  when  computing  the  intensity)  the  new  complex  amplitude 


is  found  via  the  Fresnel  transform  : 


U{io,y0)=  / / t'(xi,i/i)exp  ^r(z2+r)  hp — —-(zo-'a+JW  1)  dxxdy{  (2.15) 


Here  U{xuyx)  is  now  only  a  function  of  the  shape  of  the  exit  pupil,  in  this  case  a 


lens  (a  circle).  The  limits  of  the  integral  then  become  the  limits  of  the  aperture. 


Given  in  Figure  6  is  the  intensity  pattern  of  a  circular  aperture  containing  a  per¬ 


fect  thin  lens  in  a  plane  a  focal  length  away.  This  pattern  is  called  the  Airy  pat¬ 


tern,  named  after  G.  B.  Airy  (9). 


* 


Figure  6.  The  Airv  Pattern 


« v-.T-, 


I 

i 

I 


However,  in  taking  the  paraxial  approximation.  Goodman  excluded  any 
effects  of  aberrations  as  can  be  shown  in  Geary’s  paper  (3:52-5*1).  In  excluding  all 
of  these  aberrations,  there  would  be  no  beam  deformation  as  the  angle  of 
incidence  grows.  In  the  following  analysis,  the  binomial  approximation  will  not  be 
taken.  Therefore,  the  transfer  function  for  Goodman's  thin  lens  becomes: 


x,y)  =  exp 

r 

jk(n  — l) 

r 

[,  {*2+y2) 
R? 

1 

'  /  O  Ov  ' 

,  (*  +r) 

Rt 

1 

o 

‘) 

When  this  equation  is  used  to  find  the  field  amplitude  using  the  Fresnel 
transform,  the  intensity  pattern  is  derived  and  given  below  in  Figure  7. 


Note  that  the  main  aberration  affecting  the  system  is  spherical  aberration. 
The  second  ring  of  the  pattern  is  .0190  times  the  height  of  the  central  maximum 
whereas  in  an  unaberrated  case  the  second  ring  would  be  .0175  times  the  central 


maximum.  As  long  as  the  plane  wave  is  parallel  to  the  optic  axis,  the  contribu¬ 
tions  due  to  astigmatism  and  coma  are  zero. 

Goodman  handles  an  off-axis  plane  wave  as  a  function  of  the  input  U(x.y). 
Figure  8  below  shows  the  system  used. 


Figure  8.  The  General  Thin  Lens  Configuration 
The  input  plane  wave  then  becomes: 

Ut(x  ,y)  =  exip(jksinay)  (--17) 

versus  one  for  a  unit  amplitude  plane  wave  parallel  to  the  optic  axis,  and  alpha  is 

the  angle  in  the  y-z  plane  with  respect  to  the  z  or  optic  axis.  The  complex  ampli¬ 
tude  can  still  be  solved  with  the  methodology  given  previously.  The  complex 
amplitude  in  the  image  plane  was  solved  by  taking  the  Fast  Fourier  Transform 
(FFT)  of  the  plane  wave  incident  on  a  thin  lens.  The  output  of  the  FFT  was 
then  squared  to  obtain  the  intensity  distribution  in  the  image  plane.  The  program 
used  for  this  analysis  is  located  in  Appendix  A.  The  intensity  patterns  for  various 
angles  alpha  are  shown  in  Figure  9  for  Lens  1  and  in  Figure  10  for  Lens  2.  The 
Streh!  Ratios  for  the  intensity  patterns  decrease  as  the  angle  alpha  increases. 
Here,  the  Goodman  case  analyzed  still  carries  all  of  the  contributions  due  to 
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spherical  aberration  and  defocus,  and  as  can  be  seen  in  Chapter  4.  as  the  angle  is 
increased,  the  Zernike  coefficient  for  spherical  aberration  increases  a  finite  amount 
due  to  the  changing  defocus  term.  Therefore  the  Strehl  Ratio  should  correspond¬ 
ingly  decrease  by  a  finite  amount  as  is  shown  in  the  two  figures. 

As  is  shown  in  the  previous  graphs,  the  plane  wave  at  an  arbitrary  angle 
alpha  just  moves  the  a  =  0  pattern  up  and  down  the  image  plane  with  respect  to 
the  magnitude  of  alpha,  i.e.,  the  larger  the  angle  alpha,  the  larger  the  displaement 
in  the  image  plane  from  the  optic  axis.  Therefore  in  the  Goodman  case,  the  angle 
alpha  does  not  affect  the  amount  of  aberrations  at  all,  the  aberrations  are  all  con¬ 
tained  in  the  transfer  function  within  the  square  root.  From  this  result  there 
whould  be  no  beam  degradation  at  all  as  the  field-of-view  is  increased.  However, 
this  is  known  not  to  be  true,  based  on  experimental  evidence  obtained  using  PHA- 
SAR  at  the  Air  Force  Weapons  Laboratory  (7).  Therefore,  this  description  of  a 
thin  lens  must  be  incorrect. 
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10.  Goodman  Analysis,  Lens  2.  for  a  —  0.  5,  and  10  Degrees 
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III.  A  New  Thin  Lens  Expression  (The  Walton  Thin  Lens) 


This  chapter  explores  the  expression  of  a  thin  lens  which  is  applicable  for  off- 
axis  projections  such  as  those  used  in  multi-aperture  systems.  An  expression  is 
derived  and  then  computationally  analyzed.  Following  this  are  several 
verifications  of  the  expression's  accuracy.  The  Strehl  ratios  for  all  of  the  cases  are 
presented  as  well  as  the  limiting  case  at  a  =  0  which  should  approach  an  Airy 
pattern.  Also  included  are  the  cross-sections  at  various  angles  to  determine  which 
aberrations  are  present  and  which  aberrations  cause  most  of  the  beam  deforma¬ 
tion  as  the  angle  alpha  grows. 

m.l.  The  Mathematical  Expression 

Figure  11  describes  the  system  employed  in  the  following  analysis.  The 
center  of  the  coordinate  system  is  at  the  "center"  of  the  lens,  or  where  the  two 
halves  of  the  lens  meet  at  the  optic  axis.  The  only  approximation  made  in  this 
case  is  that  the  lens  is  a  thin  lens,  i.e.,  that  one  can  ignore  the  refraction  of  light 
at  the  lens  surface.  To  make  the  mathematics  managable,  alpha,  the  angle  of  the 
incident  plane  wave  with  the  optic  axis,  is  taken  to  be  with  respect  to  the  y-z 
plane  only.  The  thin  lens  approximation  also  states  that  throughout  the  system, 
the  plane  wave  will  stay  in  a  specific  y-z  plane  for  any  given  value  of  x. 

The  next  part  of  the  analysis  is  purely  geometrical,  with  many  of  the  details 
omitted.  First,  the  lens  thickness  is  solved  for  as  in  Chapter  2.  This  new  thick¬ 
ness  is  substituted  in  for  V  in  Eq.  (2.12)  giving  the  transmittance  function  for  the 


Figure  11.  Geometry  Used  for  Walton  Thin  Lens 


new  lens  expression. 


To  solve  for  the  distance  between  the  two  surfaces  of  the  lens,  the  equations 


of  the  two  spheres  and  the  equation  for  the  line,  which  is  two-dimensional  since 


the  plane  wave  is  only  in  the  y-z  plane,  are  equated.  To  solve  for  the  point  at 


which  the  entering  ray  exits  the  lens,  equate,  for  the: 


Where  the  point  (y,;0)  = 


Sphere  :  (V+j2)2+x| =  R$  (3.1) 

Line  :  y2  =  —  tana(z2—  z0)+y  (3.2) 

(y,0),  so  that  or  the  reference  frame  of  the  lens  becomes 


the  x-y  plane.  Also: 


Solving  these  two  equations  yields  a  quadratic  equation  in  ;2: 


z2  (sec‘a)+2z2(  —  H2  j“  —  ytan  q)+  z2  +  y2— H2  =  0 


$Zo  -  0 


where 


L 


iz  3  nsv  .<•.->  see  =5&r 


i 

I 

8 
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3  =  sec "a 


f_  =  |/?;  — //""  j"  —  ytana 
»?  =  zS+r-H2 


(3.6) 


Solving  for  the  point  (y2,z 2)  : 


(3.7) 


y  o  —  — jo  tana  -f-  y 

In  a  similar  manner,  the  following  three  equations  are  solved  for  the  point 


( it  i ,  - 1 ); 


Sphere  :  (r,— V* )2+i{  +y{  =  R{ 
Line  j  =  -  tan  q(zj— ;0)+y 


(3.8) 

(3.9) 


V'  =  (flf-//2)2  (3.10) 

Where  again,  the  reference  of  the  x-y  axis  through  the  origin  was  used  such  that 

the  point  (y.:0)  —  (y,0).  Solving  these  three  equations  as  before  gives  the  point 

(yx,z |)  with  respect  to  y  and  alpha: 


Z 1  = 


?+ 

t 

f+ 

2 

7 

u 

k 

s 

/*  . 

where: 


3<i  =  -Z[tana+y 


3  —  se<ra 

l 


(3.11) 

(3.12) 


=  |/?2  —  If-  j"  -fj/tana  (3.13) 

o  o  r  ro 

7  =  *1  +y~-H~ 

Therefore.  to  solve  for  the  distance  between  the  two  surfaces  as  a  function  of 
the  system  parameters  (Ri.Rn,  and  H)  and  the  givens  (x.  y.  and  a)  the  distance  for¬ 
mula  is  needed.  The  distance  is  then  given  by: 


2d 


diPi.P:)  =  yz-yf 


(3.1-4) 


This  form  is  two-dimensional  in  y  and  z  because  of  the  thin  lens  approximation. 
Substituting  in  for  the  two  points  derived  above,  the  distance  formula  becomes: 


_1_ 

d  =  cosq( 


(3.15) 


where 


0  2 
(3  =  sec  a 


i 

— H“  j"  —  y  tanor 
—  ff2  j“  +  ytanor 


(3.16) 


r)  =  jt+jT-H2 

Therefore  the  new  transfer  function  for  the  thin  lens  becomes: 

tt(z,y)  =  exp(jk{n-l)d)  (3.17) 

The  constant  phase  factor  was  again  ignored  since  only  intensity  patterns  are 

being  considered. 


III. 2.  Computational  Analysis 

The  intensity  patterns  were  solved  by  taking  the  FFT  of  a  plane  wave 
incident  on  a  thin  lens.  The  output  of  the  FFT  was  then  squared  to  obtain  the 
intensity  distribution  in  the  image  plane.  The  program  ALPHA  written  for  this 
purpose  is  included  in  Appendix  B.  The  thin  lens  was  described  as  a  set  of  dis¬ 
tances  from  which  a  set  of  optical  path  differences  was  defined.  In  Figure  12  are 
two  graphs  which  depict  the  distances  the  light  travels  through  the  lens.  Part  (a) 
shows  the  set  of  distances  for  an  angle  of  o  =  0:  part  (b)  the  set  of  distances  for 
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Next,  the  Strehl  Ratios  were  calculated  for  all  of  the  intensity  patterns.  As 
can  be  seen,  for  the  thin  lens  the  Strehl  ratios  correlate  well  with  information 
given  by  Mills  (10).  A  significant  problem  deduced  from  the  Strehl  Ratios  is  that 
the  power  in  the  central  peak  has  decreased  by  40  percent  as  the  angle  alpha  has 
increased  to  two  degrees  in  Lens  1.  Lens  2  shows  that  there  is  a  valid  problem 
with  the  approximation  made  when  the  lens  varies  at  all  from  being  "thin".  This 
correlates  well  since  Lens  2  was  designed  to  be  a  thick  lens.  The  distance  formula 
derived  above  is  therefore  inappropriate  for  a  lens  of  this  type. 

The  last  verification  made  was  to  vary  the  size  of  the  aperture  stop  to  show 
how  the  aberrations  change.  The  smaller  the  aperture  stop,  the  more  like  an  Airy 
pattern  the  diffraction  pattern  becomes.  The  larger  the  aperture  stop,  the  more 
aberrations  there  are  present  in  the  diffraction  pattern.  This  is  shown  in  Figure 
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Figure  15.  Aperture  Stop  =  .5.  1.0,  1.5  cm  for  Lens  1  at 


5  Degrees 


IV.  Analysis  Using  FALCON,  a  Geometrical  Ray  Trace  Program 


In  the  previous  chapter,  the  Walton  Thin  Lens  expression  was  analyzed  for 
on-axis  cases.  In  this  chapter,  the  expression  is  compared  to  the  geometrical  solu¬ 
tion  for  off-axis  projections.  The  geometrical  ray  trace  program  FALCON  was 
used  to  derive  the  aberration  coefficients  which  were  used  in  a  diffraction  analysis 
and  subsequently  compared  to  the  Walton  expression.  The  program  written  to 
compute  the  diffraction  analysis  is  included  in  Appendix  C.  The  intensity  plots 
derived  for  this  analysis  are  given  at  the  end  of  the  chapter. 


IV.l.  FALCON 


A  comparison  was  made  between  the  Walton  Thin  Lens  and  the  results  from 
the  geometrical  ray  trace  program  FALCON  (11).  Initially  the  lens  characteristics 
were  first  entered  into  FALCON,  and  since  a  geometrical  analysis  requires  an 
object,  an  object  of  variable  height  was  place  at  "infinity"  so  that  a  plane  wave  of 
various  angles  could  be  modelled.  The  results,  the  aberration  coefficients,  were 
then  computed  and  inserted  into  a  diffraction  analysis.  The  aberration  coefficients 
for  piston,  tilt,  and  defocus  were  included  in  the  diffraction  analysis  for  complete 
accuracy,  but  are  not  included  below  because  the  emphasis  in  this  thesis  is  the 
beam  deformation  due  to  the  third  order  aberrations.  Listed  below  are  the 
approximate  coefficients  for  the  two  lenses  used  in  this  analysis.  As  can  be  seen 
from  the  intensity  patterns  in  Figures  16  and  17,  the  aberrations  occured  at  a 
much  smaller  angle  than  thev  did  in  the  Walton  analvsis.  However,  the  two  ana- 
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lyses  definitely  show  similarities  in  the  shape  of  the  resultant  curves  and  in  the 
aberrations  which  dominate  as  the  angle  alpha  grows.  Table  III- 1  lists  the  various 
aberration  coefficients  for  the  two  lenses  at  various  angles  alpha. 

As  can  be  seen  from  the  intensity  plots  for  the  diffraction  analysis,  the  plots 
for  the  thin  lens  correlate  with  those  from  the  Walton  analysis.  Because  of  this 
correlation,  these  coefficients  are  a  close  approximation  to  the  actual  coefficients. 
However,  the  thick  lens  did  not  correlate  as  well  in  this  case,  so  the  coefficients  for 
the  second  lens  are  not  a  good  approximation.  This  is  due  to  the  program  being 
inadequate  in  the  limit  of  the  thick  lens.  Therefore,  the  intensity  plots  for  the 
thick  lens  will  not  be  representative  of  the  actual  intensity  distributions. 

From  the  analysis  of  FALCOX,  the  diffraction  analysis  seems  to  break  down 
as  the  angle  alpha  grows  above  ten  degrees  for  the  thin  lens,  and  five  degrees  for 
the  thick  lens.  This  is  delineated  in  the  next  chapter  as  the  aberration  coefficients 
grow  out  of  proportion  to  the  experimental  analyses. 
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Table  III- 1  Aberration  Coefficients  for  Lens  1  and  Lens  2 


« 


A* 


LENS  1 


Angle 

45  Astig. 

90  Astig. 

X  Coma 

Y  Coma 

Spherical 

0 

0.0 

0.0 

0.0 

0.0 

-0.0424 

2 

0.1744 

0.0 

0.0 

0 

0.0735 

5 

1.1094 

0.0 

0.0 

0.1931 

-0.0444 

10 

4.7075 

0.0 

0.0 

0.4541 

-0.0506 

15 

11.6450 

0.0 

0.0 

0.8712 

-0.0630 

20 

23.5473 

0.0 

0.0 

1.5881 

-0.0873 

LENS  2 


Angle 

45  Astig. 

'90  Astig. 

X  Coma 

Y  Coma 

Spherical 

0 

0.0 

0.0 

0.0 

0.0 

-0.5387 

2 

.3198 

0.0 

0.0 

0.6530 

-0.5478 

5 

2.0326 

0.0 

0.0 

1.6826 

-0.5715 

10 

8.6280 

0.0 

0.0 

3.7407 

-0.6643 

15 

21.2976 

0.0 

0.0 

6.6225 

-0.8462 

20 

42.8664 

0.0 

0.0 

11.0306 

-1.1836 

All  of  these  coefficients  have  the  units  of  waves  or  optical  path  difference. 


33 


Maximum  Intcnsi  l.q  cil  l  1  29 , 1  <rl9  ) 

Slrehl  Rulco  ~  0.9900 


Maximum  I  nt.onsl.  I  ij  uL  I]  *>*>,?.■]  ,\ 
Strehl  KuLlo  -  0.9500 


floxinuim  I  atonal  l  ^  ci L  (  20‘1  ,  1  Sf  ) 
Strehl  Ratio  -  0.1000 


mm 

dWW 

v- 


Figure  16.  FALCON  Analysis,  Lens  I ,  for  a  —  0.  5,  and  10  Degrees 


V.  Experimental  Verification 


8 

£ 

t 

l 


The  most  important  verification  for  this  new  expression  is  a  comparison  to 
experimental  results.  Below  is  the  experimental  set-up  which  was  used  to  verify 
the  impulse  response  of  a  lens.  The  lens  was  turned  through  different  angles  to 
simulate  radiation  which  would  impinge  the  static  lens  from  various  angles  in 
reality.  The  experimental  results  versus  the  computational  results  are  outlined 
below. 
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V.l.  The  Experiment 

Figure  18  gives  the  experimental  set-up  used  to  validate  the  previous 
analysis. 

The  two  lenses  described  in  Chapter  2  were  used  in  this  set-up.  As  a  rem¬ 
inder,  the  first  lens  was  a  thin  lens  with  dimensions  of  /?,  =  16.8-4  cm.  R2  =  -44.23 
cm.  H  =  4.37  cm,  and  an  aperture  stop  with  a  radius  of  1  cm.  The  second  lens 
had  dimensions  of  Rl  —  11.37  cm,  R2  —  -11.63  cm.  H  =  6.09  cm,  and  an  aperture 
stop  with  a  radius  of  1  cm.  A  HeNe  laser  was  used  in  the  experiment  and  in  all  of 
the  previous  analyses  as  well.  True  to  assumptions,  the  thicker  the  lens,  the 
worse  the  aberrations  became  as  the  angle  of  incidence  of  the  radiation  grew. 

V.2.  Results 

The  results  of  the  experiment  are  shown  in  Figures  19  and  20.  These  inten¬ 
sity  patterns  are  not  normalized  to  any  size  standard,  as  the  analysis  is  purch 


50  mW  HeNe  Laser 


Image 

Plane 


Figure  18.  Experimental  Set-Up 

qualitative  not  quantitative.  The  larger  the  off-axis  angle,  the  larger  the  pattern 
in  the  far-field.  The  distance  to  the  image  plane  is  given  for  each  set  of  pictures. 
A  correct  size  approximation  can  be  made  based  on  the  distance  from  the  tele¬ 
scope  to  the  camera,  which  was  chosen  in  order  to  get.  the  most  information  per 
picture.  All  pictures  were  taken  with  the  same  shutter  speed  and  f  number,  how¬ 
ever,  various  neutral  density  filters  were  required.  The  magnitudes  of  these  filters 
are  also  given  in  the  above  figures.  These  magnitudes  are  representative  of  the 
magnitudes  of  the  intensity  in  the  central  peak.  The  larger  the  neutral  density 
filter,  the  larger  the  intensity  value.  As  can  can  be  seen  from  the  values  of  dis¬ 
tance  and  neutral  density  filter  used,  there  is  a  much  higher  intensity  value  in  the 
central  peak  for  on-axis  radiation  than  for  off-axis  radiation.  This  is  consistent 
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with  the  Strehi  Ratio  argument  from  Chapter  3.  This  also  shows  that  one  of  the 
main  problems  with  a  large  field-of-view  is  not  the  addition  of  side  lobes,  but  the 
loss  of  power  in  the  central  peak. 

Figures  13  and  14  give  the  results  from  the  computational  analyses  of  the 
lenses  used  in  the  above  experiment.  The  data  is  normalized  to  1  cm  per  17 
matrix  elements.  As  can  be  seen,  the  plots  for  various  values  of  alpha  are  con¬ 
sistent  with  the  information  obtained  from  the  experiment. 

The  only  aberration  to  be  concerned  with  for  a  =  0  is  spherical  aberration. 
The  first  lens  was  of  poor  quality,  an  Airy  pattern  never  being  achieved  using  this 
lens.  Therefore  this  lens  exhibits  more  aberrations,  such  as  astigmatism,  than  the 
computations  would  handle.  The  experimental  results  for  the  second  lens  are  very 
similar  to  the  results  found  from  the  computational  analysis  for  a  =  0.  As  the 
incident  angle  grew  for  the  first  lens,  astigmatism  became  the  dominant  contribu¬ 
tor  to  the  intensity  pattern.  This  occured  through  an  angle  of  ten  degrees.  .After 
ten  degrees,  the  dominant  contributor  became  coma  and  at  this  point  there  is  a 
serious  problem  with  the  intensity  pattern  shifting  in  the  image  plane.  As  can  be 
expected  from  a  thicker  lens,  the  results  are  just  a  worse  case  of  the  first  lens. 
Astigmatism  became  a  problem  at  an  angle  of  two  degrees,  and  coma  became  a 
problem  at  an  angle  of  five  degrees.  As  can  be  seen  in  the  last  picture  for  the 
second  lens  (a  =  20),  the  intensity  pattern  became  so  large  that  the  telescope 
became  the  limiting  aperture.  Again,  the  results  from  the  computational  analysis 
did  not  have  this  problem.  However,  the  program’s  results  are  not  as  accurate  for 
this  second  lens  as  they  were  for  the  first.  Also,  the  matrix  size  was  too  small  for 


any  results  at  an  angle  of  alpha  greater  than  ten  degrees,  so  no  accurate  comparis¬ 
ons  can  be  made  between  analysis  and  experiment  in  this  region.  In  the  region  0 
<  a  <  10  degrees  however,  there  is  good  correlation,  just  as  had  been 
hypothesized. 

Therefore,  as  has  been  seen  from  the  experimental  results  given  above. 
Goodman's  analysis  is  a  good  approximation  for  on-axis  radiation  effects.  If 
Geary  and  Peterson's  term  is  used,  the  result  becomes  extremely  accurate  for  on- 
axis  radiation,  but,  does  not  hold  for  off-axis  radiation.  The  most  obvious 
difference  between  these  previous  results  and  the  new  expression  for  the  thin  lens 
derived  herein  is  that  instead  of  off-axis  effects  simply  moving  the  intensity  pat¬ 
tern  in  the  image  plane,  these  effects  actually  do  not  move  the  pattern  as  much  as 
cause  deformation  of  the  pattern  and  severe  aberrations  for  large  angles. 
Although  beam  deformation  does  seem  to  be  a  problem  for  the  field-of-view 
required  by  SDI,  the  amount  of  normalized  intensity,  or  power,  lost  in  the  central 
lobe  as  the  angle  alpha  is  increased  can  also  cause  serious  problems  for  detection. 


VI.  Application:  Multi- Aperture  Systems 


In  this  chapter,  muli-aperture  systems  will  be  applied  to  both  the  thin  and 
the  thick  lenses  investigated  previously,  and  the  effect  of  aberrations  on  field-of- 
view  will  be  investigated.  Note  that  in  all  of  the  multi-aperture  intensity  pat¬ 
terns,  the  envelope  of  the  pattern  is  the  intensity  pattern  for  the  single  lens  used 
in  the  system.  In  all  of  the  following  analyses,  the  optic  axis  is  located  at 
(129,129)  and  17  matrix  element  are  taken  to  be  1  cm.  The  configurations  to  be 
used  in  this  chapter  are  given  in  Figure  2.  The  program  written  to  compute  the 
intensity  patterns  for  the  multi-aperture  configurations  is  included  in  Appendix  D. 


VI, 1,  Single  versus  Multi- Aperture  Patterns 

Given  in  Figure  21  are  the  intensity  patterns  for  the  single  thin  Lens  1  at  an 
angle  a  =  5  degrees  and  the  four  aperture  case  using  this  lens  at  an  angle  of  five 
degrees  as  well.  Note  that  the  pattern  for  the  single  lens  is  the  envelope  for  the 
multi-aperture  pattern.  This  result  is  in  keeping  with  the  theory  as  shown  by 
Young  (12:281-283)  in  his  classic  double  slit  experiment. 

VI. 2.  Multi- Aperture  System  Using  Goodman  Analysis 

As  a  comparison  to  the  previous  theory,  the  Goodman  analysis  was  applied 
to  the  three  multi-aperture  systems  given  in  Figure  2.  As  can  be  seen  from  Fig¬ 
ures  22  through  2-1.  as  the  angle  alpha  increases  to  twenty  degrees  the  side  lobes 
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Figure  21.  Single  vs.  Multi-Aperture 


grow  asymmetrically  in  the  pattern.  As  the  angle  increases  through  ten  degrees 
the  side  lobes  approach  a  maximum,  and  as  the  angle  approaches  twenty  degrees 
the  side  lobes  decrease  again.  This  phenomenon  is  apparent  in  all  of  the  multi¬ 
aperture  configurations  using  the  Goodman  thin  lens  model.  The  form  of  these 
results  are  not  consistent  with  standard  results  using  Goodman’s  thin  lens.  How¬ 
ever,  the  analyses  of  interest  in  this  thesis  were  based  on  the  new  lens  equation  so 
this  discrepency  was  not  investigated.  The  only  lens  used  in  this  analysis  is  Lens 
1  since  the  patterns  for  this  lens  and  Lens  2  are  essentially  the  same.  Note  that 
the  Strehl  ratios  for  the  Goodman  analysis  decreased  by  a  finite  amount  as  in  the 
single  aperture  case.  Again,  this  is  due  to  the  effects  of  spherical  aberration  and 
defocus  as  described  in  Chapter  4. 
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VI. 3.  Three  Apertures 
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The  following  analysis  is  for  Lens  1  and  Lens  2  in  the  three  aperture 
configuration.  Given  in  Figure  25  and  Figure  26  are  the  graphs  for  the  intensity 
patterns  at  various  values  of  alpha  for  both  lenses.  .-\s  can  be  seen  from  Lens  1. 
the  central  peak  is  much  more  narrow  than  for  the  single  lens  case.  Therefore,  two 
point  resolution  for  multi-aperture  systems  is  much  better  than  for  a  single  lens  as 
has  been  reported  previously  (5:40—16).  However,  as  the  incident  radiation 
impinges  the  lens  at  an  angle  greater  than  zero,  the  side  lobes  on  one  side  of  the 
central  peak  grow  more  than  those  on  the  other  side.  This  is  apparent  from  the 
intensity  patterns  for  the  single  lens  which  form  the  envelope  for  these  plots.  At 
angles  larger  than  five  degrees  coma  is  the  dominant  aberration,  which  causes  an 
asymmetric  shift  of  the  light  intensity  to  one  side  of  the  pattern.  Thus,  there  is 
consistency  with  patterns  shown  in  Figures  13  and  14  and  the  information  given 
above.  As  the  angle  of  incidence  grows  to  five  degrees,  the  side  lobes  grow  to 
approach  the  height  of  the  central  peak.  Ii  addition,  there  is  a  small  shift  in  the  x 
axis  caused  by  astigmatism  which  has  to  be  included.  This  causes  problems  with 
detection  which  has  not  been  dealt  with  in  any  of  the  research  to  date.  However, 
with  any  kind  of  threshold  level  detector  this  problem  can  be  overcome  (5:63-61). 
From  the  graphs  in  Figure  25.  as  alpha  grows,  the  main  lobe  separates  into  many 
peaks,  causing  definite  problems  for  two  point  resolution  as  pointed  out  by  Wat¬ 
son  (5:10-62).  An  additional  problem  for  this  optical  system  is  that  the  amount 
cf  power  in  the  central  peak  at  an  angle  of  five  degrees  is  approximately  10  per¬ 
cent  < f  that  in  the  o  —  0  case. 
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Lens  2,  as  seen  in  Figure  26.  is  a  poor  lens  based  on  multi-aperture  considera¬ 
tions.  At  an  angle  of  zero  degrees,  the  side  lobes  of  the  system  are  approximately 
45  percent  of  the  intensity  of  the  central  peak.  This  is  obviously  not  a  good 
multi-aperture  system.  As  the  angle  grows,  the  side  lobes  increase  to  the  point 
that  by  the  time  the  angle  is  two  degrees,  the  side  lobes  are  within  90  percent  of 
the  central  peak.  At  an  angle  of  ten  degrees  the  central  intensity  begins  to  walk 
down  the  x  as  well  as  the  y  axis.  Again,  this  must  be  accounted  for  in  the  detec¬ 
tion  system.  For  angles  above  two  to  five  degrees,  virtually  no  information  can  be 
gleaned  from  this  system  with  the  thick  lens  because  of  the  many  peaks  in  the 
central  region.  In  addition,  the  central  intensity  has  decreased  to  approximately 
one  percent  of  its  initial  value.  The  aberrations  present  as  the  angle  increases 
indicate  that  this  lens  is  unacceptable  for  multi-aperture  systems. 

VT.4.  Four  Apertures 

In  Figure  27  are  the  intensity  patterns  for  the  four  aperture  system  at  vari¬ 
ous  angles  alpha.  The  four  apertures  almost  completely  fill  a  square  in  the  optical 
system,  thus  a  square  pattern  is  achieved  in  the  image  plane.  The  multi-aperture 
intensity  plots  for  Lens  1  are  much  better  than  for  Lens  2  based  on  the  results 
from  section  YI.3.  Therefore,  the  analysis  of  Lens  2  for  the  four  aperture  case  will 
be  omitted.  For  Lens  1,  the  side  lobes  are  very  small  at  n  =  0.  However,  as  the 
angle  approaches  two  degrees,  the  side  lobes  are  approximately  85  percent  of  the 
central  peak  intensity.  Again,  as  the  angle  alpha  increases,  the  number  of  peaks  in 
the  same  range  as  the  central  maximum  increase,  and  thus  there  are  problems 
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Figure  25.  Walton  Lens  1  in  Three  Aperture  Configuration 
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with  a  large  field-of-view  for  this  system.  As  before,  the  central  peak  shifts  down 
the  x  axis,  by  a  relatively  small  amount  compared  to  the  shift  down  the  y  axis. 

VI. 5.  Six  Apertures 

Figure  28  gives  the  intensity  patterns  for  the  six  aperture  system  at  various 
angles  alpha.  As  can  be  seen  from  the  patterns  for  Lens  1.  the  central  peak  has 
decreased  in  width  significantly  from  the  three  aperture  case,  and  the  side  lobes 
are  also  much  smaller  than  in  any  of  the  previous  systems.  These  are  excellent 
results  based  on  two  point  resolution  criteria.  However,  as  the  angle  alpha  grows, 
many  side  lobes  increase  to  a  height  equivalent  to  the  central  peak  by  an  angle  of 
five  degrees.  This  optical  arrangement  actually  proves  to  be  a  far  worse 
configuration  than  that  for  the  three  aperture  case.  This  result  contradicts  previ¬ 
ous  research  (5,6).  Therefore,  this  optical  system  shows  the  same  problems  as  the 
previous  two  optical  systems.  This  is  an  unexpected  result  based  on  previous 
research  concerning  SDI’s  required  twenty  degree  field-of-view.  No  comment  will 
be  made  about  the  six  aperture  case  for  Lens  2  given  the  results  from  the  three 
aperture  case. 

VI. 6.  Field-of-View  and  Power  Considerations 

Field-of-view  and  power  in  the  central  peak  are  the  two  basic  problems  of 
interest  within  this  thesis.  As  can  be  seen  from  Figures  25  through  28,  the  field- 
of-view  of  any  system  varies  with  the  number  of  apertures  in  the  system,  their 
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placement,  and  the  angle  of  Incidence  cf  the  incoming  radiation.  Only  Lens  1  will 


be  considered  in  this  synopsis  since  the  intensity  patterns  for  Lens  2  were  so  poor 
when  compared  to  Lens  1.  As  determined  previously,  the  more  apertures  there  are 
in  the  optical  system,  the  narrower  the  central  peak,  which  i.s  a  very  good  result 
for  two  point  resolution.  Also  aperture  placement  is  important  in  field-of-view 
considerations.  While  the  four  aperture  case  was  poor  when  compared  to  the 
three  aperture  case,  the  six  aperture  case  seemed  to  be  the  best  of  all  systems  stu¬ 
died.  However,  this  is  apparently  only  for  on-axis  radiation.  The  computations 
achieved  in  this  thesis  indicate  that  for  olT-axis  radiation  with  a  more  accurate 
analysis  for  the  transfer  function  of  a  thin  lens,  the  beam  degradation  actually 
becomes  far  worse  as  the  number  of  apertures  is  increased.  Also,  there  exists  a 
problem  of  threshold  levels  dependent  upon  the  angle  of  incidence.  At  angles  less 
than  five  degrees  these  systems  behave  well,  but  for  angles  greater  than  five  to  ten 
degrees  the  system  performance  degrades  due  to  the  large  number  of  side  lobes  of 
the  same  approximate  size  as  the  central  peak.  The  structure  of  these  patterns, 
i.e.  the  envelope,  is  the  same  as  that  obtained  in  the  experiment. 

The  main  problem  forseen  in  future  optical  systems  is  not  the  problem  of 
field-  of-view  in  relation  to  two  point  resolution  but  the  problem  that  the  further 
off-axis  the  incident  radiation,  the  le:  power  appearing  in  the  central  lobe.  This 
could  cause  serious  problems  with  tection  as  has  been  mentioned  previously.  If 
insufficient  intensity  exists  to  detect  the  signal,  aberrations  are  essentially  unin- 
portant.  Also,  if  any  kind  of  threshold  system  were  used  for  this  detector,  the 
angle  of  incidence  must  be  considered  so  that  the  threshold  could  be  varied.  This 


poses  an  extremely  difficult  design  problem  for  future  optical  systems. 
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VII.  Results:  A  Comparison  of  the  Four  Methods 


Given  in  this  chapter  is  a  comparison  of  the  four  analyses  discussed  previ¬ 


ously  to  find  an  overall  method  for  describing  a  transfer  function  of  a  th 


in  lens 


appropriate  for  off-axis  projections.  The  four  methods  discussed  are:  Goodman’s 


Analysis.  The  Walton  Analysis.  FALCON,  and  Experiment. 


VII.l.  Comparison  of  the  Results 


Given  in  Figures  29  through  32  are  a  comparison  of  the  results  at  the  specific 


angles  zero.  two.  five,  and  ten  degrees  for  Lens  1.  At  zero  degrees,  all  four  of  the 


methods  correlate  extremely  well.  This  is  the  limiting  condition  for  the  Walton 


analysis,  and  the  validity  regime  for  the  Goodman  analysis.  This  result  is  excel¬ 


lent  since  it  shows  consistency  with  the  results  shown  previously  in  the  literature 


(1,8).  The  only  aberration  present  in  this  analysis  is  spherical  aberration. 


At  an  angle  of  two  degrees,  the  Goodman  analysis  shows  that  there  is  no 


beam  deformation,  only  a  translation  in  the  y  axis  from  the  plane  wave.  However. 


the  experiment  explicitly  shows  how  astigmatism  is  present  as  well  as  a  smal 


translation  in  the  y  axis.  FALCON  shows  not  only  an  astigmatic  problem,  but 


also  a  contribution  due  to  coma.  The  Walton  analysis  shows  the  same  intensity 


listribution  as  the  experiment. 


At  an  angle  of  five  degrees,  the  Goodman  analysis  again  shows  only  a  trans¬ 


lation  in  the  y  axis.  The  experiment  shows  not  only  an  astigmatic  problem,  but  a 


large  contribution  due  to  coma.  The  intensity  distributions  from  FALC'OX  show 
only  the  portion  of  the  coma  pattern  in  the  v  direction,  the  point  of  maximum 
intensity.  The  Walton  analysis  at  this  point  .hows  not  only  the  portion  of  the 
coma  pattern  in  the  v  direction,  but  the  astigmatic  asymmetry  in  the  x  direction 
as  well.  The  Walton  analysis  correlates  well  with  the  experimental  results. 

At  an  angle  of  ten  degrees,  the  Goodman  analysis  shows  a  further  transla¬ 
tion.  The  experimental  result  shows  a  large  component  due  to  coma,  which 
literally  washes  out  the  comtribution  due  to  astigmatism.  The  FALCON  analysis 
has  exceeded  the  limitations  of  the  program  because  of  the  finite  size  of  the  matrix 
used  in  the  FFT.  The  Walton  analysis  again  correlates  well  with  the  experimental 
results. 

Given  in  Figures  33  through  36  are  the  results  for  Lens  2.  The  analysis  is  the 
same  as  that  of  Lens  1  except  that  the  strength  of  the  aberrations  occur  at  much 
smaller  angles  than  for  Lens  1.  Also,  the  Walton  analysis  has  broken  down  since 
a  thick  lens  does  not  agree  with  the  basic  assumption  of  the  analysis.  FALCON, 
as  before,  breaks  down  after  five  degrees  since  the  analysis  has  exceeded  the  limi¬ 


tations  of  the  program. 
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Figure  30.  Lens  1  Comparison  at  a  =  2  Degrees 
(a)  Experiment,  (b)  Goodman,  (e)  Walton,  (d)  FALCON 


Maximum  InlonscLij  ot.  (  1  .50 , 1  ) 

Slrehl  Ratio  ’  ().lJ(Jri/ 


Figure  31.  Lens  1  Comparison  at  a  —  5  Degrees 
(a)  Experiment,  (b)  Goodman,  (c)  Walton,  (d)  FALCON 
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Figure  32.  Lens  1  Comparison  at  q  =  10  Degrees 
(a)  Experiment,  (b)  Goodman,  (c)  Walton,  (d)  FALCON 
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Figure  3o.  Lens  2  Comparison  at  a  =  5  Degrees 
(a)  Fxperiment.  (h)  Goodman,  (c)  Walton,  (d)  l  ALC  ON 


VIII.  Conclusions 
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This  thesis  contains  a  number  of  breakthroughs,  the  most  important  of 
which  is  a  new  expression  for  a  thin  lens  applicable  for  off-axis  projections.  This 
expression  qualitatively  analyzes  the  aberrations  induced  as  the  incident  radiation 
impinges  a  lens  at  an  arbitrary  angle  alpha  with  respect  to  the  optic  axis.  This 
expression,  unlike  Goodman's  (1:S0),  displays  how  the  resulting  aberrations  will 
affect  the  field-of-view  of  a  system.  This  system  was  verified  through  various 
methods,  the  most  important  of  which  was  a  simple  experiment.  This  experiment 
explicitly  shows  how  the  beam  deformation  in  the  image  plane  affects  the  field-of- 
view  of  an  optical  system.  Through  this  experiment,  the  accurate  modelling  of 
the  system  by  a  computer  program  was  determined. 

Subsequently,  the  program  ALPHA  was  used  to  model  two  different  lenses, 
one  thin  and  one  thick.  The  thick  lens  was  employed  to  test  the  limitations  of 
the  underlying  approximations.  This  lens  accurately  verified  the  extent  of  the 
thin  lens  approximation.  As  the  angle  alpha  grew,  the  finite  limitation  of  the 
matrix  size  in  the  program  limited  the  field-of-view  of  the  results  derived  from  the 
W  alton  analysis. 

The  thin  lens  behaved  very  well,  both  in  the  experimental  analysis,  and 
under  the  computational  analysis.  This  lens  showed  accurately  both  the  problems 
and  re-uhs  of  increasing  the  angle  alpha.  As  the  angle  was  increased  from  zero  to 
five  degrees,  astigmatism  became  the  dominant  aberration.  As  the  angle  was  then 
increased  from  five  to  ten  degrees,  astigmatism  was  still  present,  but  coma  became 
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the  dominant  aberration.  Both  of  these  aberrations  are  off-axis  aberrations  as 
defined  by  Born  and  Wolf  (8:211-218).  At  an  angle  of  zero  degrees,  the  only  aber¬ 
ration  present  was  spherical.  The  amount  of  spherical  aberration  present  stayed 
basically  constant  as  the  angle  was  increased. 

The  most  important  result  of  this  thesis  is  the  application  of  this  lens  expres¬ 
sion  to  the  multi-aperture  configurations  listed  in  Chapter  VII.  These  results  are 
very  important  for  SDI  initiatives.  A  tield-of-view  of  twenty  degrees  is  possible 
from  any  multi-  aperture  system  made  up  of  thin  lenses.  However,  the  number  of 
lenses,  placement  of  the  lenses,  and  the  angle  of  incidence  of  the  incoming  radia¬ 
tion  are  vitally  important  criteria  when  designing  the  most  accurate  optical  sys¬ 
tem.  The  more  apertures  the  better  the  field-of-view  for  on- axis  radiation.  How¬ 
ever.  for  off-axis  radiation,  as  alpha  increases,  the  fewer  apertures  the  better  for 
field-of-view  considerations.  By  comparing  the  three  and  four  aperture  case,  the 
placement  of  the  lenses  is  seen  to  be  a  dominant  factor  in  determining  the 
optimum  multi-aperture  design.  Conversely,  in  all  three  systems  modelled,  there 
was  a  problem  in  the  image  plane  of  the  side  lobes  approaching  the  height  of  the 
central  maximum  at  angles  of  incidence  between  eight  and  twelve  degrees  for  the 
(loodman  analysis.  Further  research  should  be  done  in  this  area  to  ascertain  the 
problems  inherent  in  this  angular  range  of  incident  radiation  based  on  this  specific 
type  ( ff  analy-i'. 

The  ni'M  important  result  of  this  thesis  is  that  while  the  system  fie|d-of-viev> 
is  a  siyiiitb’aiit  problem,  tlm  power  lost  in  the  central  maxin  am  is  the  largest  obs¬ 
tacle  to  overcome.  As  the  incident  angle  approached  ten  degrees,  the  power  in  the 


central  peak  decreased  to  one  percent  of  the  power  in  the  on-axi.s  case. 

A  possible  solution  to  this  aberrational  problem  lies  in  the  engineering  of  the 
overall  optical  system.  If  a  beam  steering  optic,  such  as  a  grating,  were  placed  in 
the  optical  path  of  the  incident  radiation  before  the  first  optical  element  of  the 
receiver  it  could  possibly  alter  the  path  from  off-axis  to  on-axis.  This  would 
enable  the  system  to  use  the  best  imaging  criteria^  of  a  multi-aperture  system 
while  keeping  the  advantages  of  an  optical  beam  steering  device.  More  research 
should  also  be  conducted  in  the  arenas  of:  the  application  of  incoherent  light  to 
the  new  expression  for  the  thin  lens,  the  effects  of  off-axis  projection  within  the 
angles  of  eight  to  twelve  degrees  in  the  Goodman  multi-aperture  analysis,  and  the 
design  of  accurate  multi-  aperture  systems  incorporating  this  new  lens  expression. 


APPENDIX  A 


This  appendix  contains  the  code  for  the  program  GOODMAN.  This  program 
calculates  the  far-held  intensity  pattern  of  the  impulse  response  for  the  Goodman 
thin  lens  at  various  angles  of  incident  radiation  (with  respect  to  the  optic  axis). 
The  aperture  is  defined  by  the  knowns:  the  two  radii  of  curvature  for  the  lens, 
the  radius  of  the  aperture  stop,  and  the  angle  of  incidence  a.  The  various  dis¬ 
tances  through  the  lens  are  then  calculated  for  various  x  and  y  coordinates  in  the 
central  plane  of  the  lens  as  is  defined  in  Figure  4.  These  distances  are  then  com¬ 
puted  as  phase  lags  for  the  incident  plane  wave  (at  an  arbitrary  angle).  This 
aperture  is  then  Fresnel  propagated  to  the  image  plane  where  the  intensity  plots 
are  computed  and  given  in  the  text.  The  parameters,  as  well  as  other  dimensions 
determined  by  the  computer,  are  in  units  of  pixels.  The  optic  axis  is  centered  on 
the  position  (12S.5.12S.5)  in  the  object  plane,  and  in  the  image  plane  the  system 
rounds  the  optic  axis  to  be  at  the  position  (129.129). 


I 


c 

c 

C  PROGRAM:  GOODMA.VF 

C 

C 

C  THIS  PROGRAM  CALCULATES  TILE  INTENSITY  OF  THE  TWO 
C  DIMENSIONAL  DL-YGE  OF  A  PLANE  WAT  THROUGH  A 
C  THIN  LENS  AT  SOME  ANGLE  ALPHA 
C 
C 

INTEGER  L\  1 .  L\2 .  N 1 , N2 .  N3 , 1 J OB  .1  WK(2000 ),  N 
INTEGER  LARG.J.VRG. CENTER 
RE.AL  RWK(2000).LNT(256.256).X.Y,M(25G.256).Rl.R2 
REAL  INDEX.PI.RAD.D  ARG2.F.  ALPHA. \ LAX.  ASTOP 
COMPLEX  CWK(256).A(256.256).B(256.25G) 

C 


OPEN(l'NIT=23,STATUS'='NEW’.FILE='.APER.DAT') 

OPEN(UNIT=2-l.STATUS='NE\V\FILE='DBGFFT.DAT’) 

OPEN(l’NIT=2o.STATL'S='NEW',FILE=’TEST.DAT') 

OPEN(UNIT=26,STATL*S='NEW',FILE='PERFECTA.DAT') 

C 

C  READ  PARAMETERS 


N=25G 
PI=3.1  115927 
WRITE(5.29) 

READ(5.30)  Rl 
READ(5.30)  R2 
READ(5.30)  ASTOP 
READ(5.30)  INDEX 
READ(5.30)  ALPHA 
READ(5.31)  CENTER 

FORXL\T(2X.‘TYPE  INVALIDS  FOR  Rl.R2.ASTOIMNDEX.ALiM  I  A.  AND 


C'ENTER(0  OR  1).  RETURN  .AFTER  EACH  ENTRY.  REMEMBER 

50  FORMAT(F  15.-1) 

51  FORMAT(  13) 

WRITE (5. 30)  Rl 


MODS!' 


Rl  — RADIUS  OF  CURVATURE  OF  LEFT  HAND  SURF  \<  I  . 
R2= RADIUS  OF  CURVATURE  OF  RIGHT  HAND  M  Cl  \« 
ASTOP=RADIUS  OF  APERTURE  STOP  OF  THE  sYM'l  M 
I ND E X = I NI ) E X  OF  REFRACTION  OF  THE  LEX"  Ml  DM  \ 
F— FOCAL  LENGTH  OF  THE  THIN  El  Ns 
ALPHA = ANGLE  OF  INCIDENCE  FOR  PLANE  W  \\  i 
CENTER  =  WHETHER  OR  NOT  YOU  WANT  Till  INTI  N- 
CENTERED  ON  THE  MAXIM'  XI  OR  NOT 


FIX  <  i i :  \i’i i 


a 


$ 

s| 

*» 

li 


vJi 

i 

•'» 


$ 

•V 

•s 


I 


#*CI 


* 


AL  P  HA= ALP  HA  *  P I  / 1  SO. 

F=1./((INDEX-1.)*(1/R1-1/R2)) 

DELX=AST0P/17. 

CREATE  .APERTURE 

DO  220  1=1, N 
XI=I 

DO  200  J=1,N 
YJ=J 

RAD=DELX*SQRT((XI-(l27.+0.5))**2+(YJ-(127.+0.o))**2) 

X=(XI-(127.+0.5))*DELX 

Y=(YJ-(127.+0.5))*DELX 

IF(R.AD.LT.ASTOP)  GO  TO  60 

EF(R.AI).GT.ASTOP)  GO  TO  100 

D=-Rl*(l.-SQRT(1.-((X**2+Y**2)/Rl**2)))+ 
g>  R2*(1.-SQRT(1.-((X**2+Y**2)/R2**2))) 

ARG1=PI*(X**2+Y**2)/F 

ARG2=2.*PI*(IXDEN-1.)*D+ARG1-2.*PI*X*SIX(ALPFLA.) 
IF(I.E0.128)  THEN 
\VRITE(24,87)  J.X,  Y.D ,  AR G 1  ,.AP G2 
F0RMAT(2X,I3,5(3X,1PE11.4)) 

END  IF 

A(U)=C\IPLX(COS(ARG2),SIN(XRG2)) 

GO  TO  120 
A(I.J)=(0. 0,0.0) 

CONTINUE 

CONTINUE 

CONTINUE 

OUTPUT  .APERTURE  MATRIX 

DO  230  1=90.127 
DO  225  J=90.127 

\VRITE(23.40)  I.J.A(I.J) 

FORMAT!  2X.I3.3X.I3.2(3XM  PE  11.1)) 

CONTINUE 
CONTINI  E 


PERFORM  FAST  FOURIER  TRANSFORM 


I A  1  =  N 
IA2=N 
N 1  =  N 
N2  =  N 


INC;  IMS! 


c 

c 

c 


N3=l 

IJ0B=1 

CALL  FFT3D  (A,L\l.L\2,Nl,N2.N3.IJOB.I\\TERWK,CWK) 


OUTPUT  NEW  \L4TREv  A 


DO  242  J=1.N 
1=127 

WRITE(24,40)  I.M(I,J) 

-10  FORMAT(2X,I3,3X.I3.2(3X,lPEll.4)) 

242  CONTINUE 
C 

c 

IF(CENTER.EO.l)  GO  TO  707 
IF(CENTER.EQ.O)  THEN 
C 

C  SHIFT  ZERO  FREQUENCY  TO  ARRAY  CENTER 

C 

C 

DO  260  1=1.236 
DO  255  J=  1.128 
B(I,J)=A(I.J+128) 

B(I,J+128)=A(I,J) 

255  CONTINUE 
260  CONTINUE 
DO  270  J=  1.256 
DO  265  1=1,128 
A(I.J)=B(I+128.J) 

A(I+128.J)=B(I.J) 

265  CONTINUE 
270  CONTINUE 
C 
C 

C  REDUCE  SIZE  OF  ARRAY 

C 

C 

DO  250  1=1,96 
DO  2-15  J=l,96 
B(I,J)=P  fI+80,J+80) 


I 


i  .»  *  ♦ 

£ 


215  CONTINUE 
250  CONTINUE 
C 


C 

C 

c: 


WRITE  CENTERED  A  YL-VTRIX 


DO  274.1=1,96 
DO  272  .1  =  1.96 


\YRITE(2o,40)  I.J.B(I.J) 
272  CONTINUE 
274  CONTINUE 
C 


C  CALCULATE  THE  MODULUS  AND  INTENSITY  OF  .ARRAY 

C 

C 

DO  400  1=1,96 
DO  380  J=1.96 
M(I.J)=CABS(B(I,J)) 

INT(I,J)=M(I,J)**2 
380  CONTINUE 
400  CONTINUE 
C 
C 

C  FIND  THE  POINT  OF  MAXIMUM  INTENSITY 
C 

NLAX=0.0 
DO  343  1=1,96 
DO  342  J=l,96 
IF(INT(I,J).GT.MAX)  THEN 
IMAX=I 
JNLAX=J 
MAX=INT(I,J) 

ENDIF 

342  CONTINUE 

343  CONTINUE 
C 

C 

C  OUTPUT  OF  INTENSITY 
C 

IMAX= IN  LAX+80 
JMAX=JNlAX+80 
WRITE(26,134)  L\L\X,JMAX,MAX 
134  F0RMAT(2X,I3,3X,I3,3X,lPEll.4) 

DO  404  1=1,96 
DO  402  J=  1 ,96 
WRITE(26,41)  I,J,INT(I,J) 

41  F0RMAT(2X,I3,3X,I3,3X,1PE11.4) 

4U2  CONTINUE 
404  CONTINUE 
C 
C 

ENDIF 
GO  TO  999 

£  ************************************************************ 
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c 

C  TO  GR.\PH  MAXIMUM  INTENSITY  AT  CENTER  OF  PLOT 
C 

C  SHIFT  ZERO  FREQUENCY  TO  .ARRAY  CENTER 
C 

767  DO  654  1=1,256 
DO  653  J=l,128 
B(I,J)=A(I,J+128) 

B(I,J+128)=A(I.J) 

653  CONTINUE 

654  CONTINUE 
DO  652  J=  1,256 

DO  651  1=1,128 
A(I,J)=B(I+125,J) 

A(I+128,J)=B(I,J) 

651  CONTINUE 

652  CONTINUE 
C 

C 

C  CALCULATE  THE  MODULUS  AND  INTENSITY  OF  THE  ARRAY 

C 

C 

DO  721  1=1,256 
DO  720  J=  1.256 
M(I,J)=CABS(A(I,J)) 

INT(I,J)=M(I,J)**2 

720  CONTINUE 

721  CONTINUE 
C 

C 

C  FIND  THE  POINT  OF  MAXIMUM  INTENSITY 
C 

MAX=0.0 
DO  124  1=1,256 
DO  123  J= 1,256 
IF(INT(I,J).GT.MAX)  THEN 
IMAX=I 
JMAX=J 
MAX=INT(I,J) 

ENDIF 

123  CONTINUE 

124  CONTINUE 
C 

C 

C  SELECT  BOX  TO  GRAPH 
C 

IOFFSET=IMAX-4S-l 


75 


JOFFSET=JMAX-48-l 

C 

IF(IOFFSET.LT.l)  THEN 
IOFFSET=0 
END  IF 

IF(IOFFSET+96.GT.256)  THEN 
IOFFSET=256-96-l 
END  IF 
C 

IF(JOFFSET.LT.l)  THEN 
JOFFSET=0 
END  IF 

IF(JOFFSET+96.GT.256)  THEN 
JOFFSET=256-96-l 
END  IF 
C 

DO  542  1=1,96 
DO  541  J=l,96 

M(I,J)=INT(IOFFSET+I,JOFFSET+J) 

541  CONTINUE 

542  CONTINUE 
C 

C 

C  OUTPUT  OF  INTENSITY 
C 

WRITE(26,34)  IMAX,JMAX,MAX 
34  FORMAT(2X,I3,3X,I3,3X,lPEll.4) 

DO  702  1=1,96 
DO  701  J=l,96 
\VRITE(26,73)  I,J,M(I,J) 

73  F0RMAT(2X,I3,3X,I3,3X,1PE11.4) 

701  CONTINUE 

702  CONTINUE 
999  STOP 

END 


APPENDIX  B 


This  appendix  contains  the  code  for  the  program  ALPHA.  This  program  cal¬ 
culates  the  far-held  intensity  pattern  of  the  impulse  response  for  the  Walton  thin 
lens  at  various  angles  of  incident  radiation  (with  respect  to  the  optic  axis).  The 
aperture  is  defined  by  the  knowns:  the  two  radii  of  curvature  for  the  lens,  the 
height  and  aperture  stop  for  the  lens,  and  the  angle  of  incidence  a.  The  various 
distances  through  the  lens  are  then  calculated  for  various  x  and  v  coordinates  in 
the  central  plane  of  the  lens  as  is  defined  in  Figure  11.  These  distances  are  then 
computed  as  phase  lags  for  the  incident  plane  wave  (at  an  arbitrary  angle).  This 
aperture  is  then  Fresnel  propagated  to  the  image  plane  where  the  intensity  plots 
are  computed,  and  given  in  the  text.  The  parameters,  as  well  as  other  dimensions 
determined  by  the  computer,  are  in  the  units  of  pixels.  The  optic  axis  is  centered 
on  the  position  (128.5,128.5)  in  the  object  plane,  and  in  the  image  plane  the  sys¬ 
tem  rounds  the  optic  axis  to  be  at  the  position  (129,129). 


C  PROGRAM:  ALPHA.F 

C  THIS  PROGRAM  CALCULATES  THE  INTENSITY  OF  THE  TWO 
C  DIMENSIONAL  LMAGE  OF  A  PLANE  WAVE  THROUGH  A 
C  THIN  LENS  AT  SOME  ANGLE  ALPHA 
C 
C 

INTEGER  LA1,LA2.N1.X2,N3.IJOB.IWK(2000).N 
INTEGER  F1.F2, CENTER 

REAL  RWK(2000).INT(256.256).X,Y,M(2o6.256).Rl,R2.H.ALPRA 
REAL  INDEX.PI,RAD,BETA.ETA.XSIM.XSIP.ARGl.ARG2.F.\LAX 
REAL  L. AMB D A . D  ( 2 5 6 , 2 5 6 ) . XR G 3 .  AR  G4 
COMPLEX  CWK(256),A(256.256).B(256.256) 

C 

OPEN(UNIT=23,STATUS='NEW\FILE='SURFACE.DAT') 
OPEN(UNIT=24,STATUS=’XEW’.FILE=‘DBGFFT.DAT') 
OPEN(UNIT=25,STATUS=’NEW’,FILE=LVLPIL\.DAT’) 
0PEN(UN1T=26,STATUS=’NEW\FILE=LALPRA.DAT') 
OPEN(UNIT=27,STATUS=’NEW’,FILE= AIILLSA.DAT') 

C 

C  READ  PARAMETERS 
C 

N=256 
PI=3. 1415927 
\VRITE(5,29) 

READ(5,30)  Rl 
READ(5,30)  R2 
READ(5,30)  H 
RE  AD  (5, 30)  ASTOP 
READ(5,30)  ALPRA 
READ(5,30)  INDEX 
READ(5,31)  CENTER 

29  F ORMAT(2X, ’ TIT E  IN  VALUES  FOR  Rl,R2,H, ASTOP, ALPRA. INDEX 

*  AND  CENTER,  RETURN  .AFTER  EACH  ENTRY,  REMEMBER  PERIODS!') 
31  FORM\T(I3) 

30  F0RMXT(F15.4) 

C 

C  Rl=RADIUS  OF  CURVATURE  OF  LEFT  HAND  SURFACE 
C  R2= RADII'S  OF  CURVATURE  OF  RIGHT  RAND  SURF  ACT 
C  H=HEIGHT  OF  LENS  FROM  OPTIC  AXIS  AT  CENTER  OF  LENS 
C  ALPRA = ANGLE  WITH  RESPECT  TO  OPTIC  AXIS  OF  INCIDENT 

C  PLANE  WAVE 

C  INDEX=INDEX  OF  REFRACTION  OF  THE  LENS  MEDIUM 

C  F=FOCAL  LENGTH  OF  THE  THIN  LENS 

C  LAMBDA= WAVELENGTH  USED  IN  EXPERIMENT  (HENE  =  .632SE-6) 

C 

/ALPHA = ALP  II A  *  PI  / 1 80. 

F=l. /((INDEX- l.)*(l/Rl-l/R2)) 


DELX= ASTOP  / 1 7. 

LAMBD  A= .6328 
WRITE(5,30)  ALPRA 
C 

C  CREATE  .APERTURE 
C 

DO  220  1=1, N 
XI=I 

DO  200  J=1,N 
YJ=J 

RAD=DERMSQRT((XI-(l27.+0.5))**2.+(YJ-fl27.+0.5))**2.) 
X=(XI-(127.+0.5))*DELX 
Y=(YJ-(127.+0.5))*DELX 
rF(R-AD.LT.ASTOP)  GO  TO  60 
IF (RAD . G T . AS T OP )  GO  TO  100 
60  BETA=  1./ (C  OS(A_LPHA)*  *2.) 

ETA=X**2.+Y**2.-H**2. 

XSIM=(SQRT(R2**2.-H**2.))-Y*TAN(ALPRA.) 

XSIP=(SQRT(Rl**2.-H**2.))-fY*TAN(ALPR\) 

D(I,J)=COS(-\LPHA)*((SQRT(XSIM**2.-ETA*BETA))+ 

*  (SQRT(XSIP**2.-ETA*BETA))-XSIM-XSIP) 

ARGl=PI*(X**2.+Y**2.)*10**4/(F*LAMBDA) 
ARG2=2.*PI*Y*SIN(ALPRA)*10**4/LAMBDA 
ARG3=2.*PI*(INDEX-1.)*D(I,J)*10**4/LAMBDA 
ARG4=ARG3+ARGl-ARG2 
WRITE  (27, 97)  I,J,D(I,J).ARG1.ARG2.ARG3.ARG4 

97  FORMAT(2X,2(I3,3X),5(lPEll.4,3X)) 

IF(I.EQ.127)  THEN 

WRITE(27,98)  J, RAD  ,X,Y, BETA, ETA, XSIM,XSIP 
WRITE(27,99)  D,ARGl,ARG2 

98  F0RMA.T(2X,I3,7(2X,1PE11.4)) 

99  F  ORMAT  ( 1 5X,3(  1PE 1 1 .4 .2X)) 

END  IF 

A(I,J)=CMPLX(COS(.ARG4),SIN(-ARG4)) 

GO  TO  120 

100  A(I,J)=(0. 0,0.0) 

D(I,J)=0.0 

120  CONTINUE 

200  CONTINUE 
220  CONTINUE 
C 

C  OUTPUT  -APERTURE  MATRIX 
C 

DO  230  1=104,151 
DO  225  J=  104,1 51 

IF(REAL(A(I,J)).EQ.O.)  THEN 
M(I,J)=0.0 


ELSE 

M(LJ)=ATAN(AINL\G(A(LJ))/REAL(A(I,J))) 

EXDIF 

225  CONTI  NIT 

230  CONTINUE 
C 
C 

C  PERFORM  FAST  FOURIER  TRANSFORM  USING  EMSL 
C' 

C 

210  L\1=N 

L\2=N 
N1=N 
N2=N 
N3=l 
IJOB=l 

CALL  FFT3D  (A,LUl,LC2,NT,N2,N3,IJOB,IUK,R\Mv.CW7<) 
C 

C  OUTPUT  NEW  NLUTREX  A 
C 

DO  242  J=1,N 
1=127 

WRITE(24,40)  I,J,A(I,J) 

40  FORMAT(2X.I3,3X,I3,2(3X,lPEll  .4)) 

242  CONTINUE 

C 

C 

IF(CENTER.EQ.l)  GO  TO  767 
IF(CENTER.EQ.O)  THEN 
C 
C 

C  SHIFT  ZERO  FREQUENCY  TO  ARRAY  CENTER 

C 

C 

DO  260  1=1,256 
DO  255  J=l,128 
B(I,J)=A(I,J+128) 

B(I,J+128)=A(I,J) 

255  CONTINUTE 

260  continut: 

DO  270  J= 1,256 
DO  265  1=1,128 
A(I,J)=B(I+128,J) 

A(I+128,J)=B(I,J) 

265  CONTINUE 
270  CONTINUE 
C 


c 

C  REDUCE  SIZE  OF  ARRAY 

C 

C 

DO  250  1=1.96 
DO  245  J=l,96 
B(I,J)=A(I+80.J+80) 

245  CONTINUE 
250  CONTINUE 
C 
C 

C  WRITE  CENTERED  A  MATRIX 
C 

DO  274  1=1,96 
DO  272  J=l,96 
WRITE(25,40)  I,J,B(I,J) 

272  CONTINUE 
274  CONTINUE 
C 
C 

C  CALCULATE  THE  MODULUS  AND  INTENSITY  OF  ARRAY 
C 

c 

DO  400  1=1,96 
DO  380  J=  1 ,96 
M(I,J)=CABS(B(I,J)) 

INT(I,J)=M(I.J)**2 
380  CONTINUE 
400  CONTINUE 
C 
C 

C  FIND  THE  POINT  OF  MAXIMUM  INTENSITY 
C 

MAN=0.0 
DO  343  1=1,96 
DO  342  J=  1,96 
IF(INT(I,J).GT.MAX)  THEN 
L\LAX=I 
JMAX=J 
MAX=INT(I,J) 

ENDIF 

342  CONTINUE 

343  CONTINUE 
C 

C 

C  OUTPUT  OF  INTENSITY 
C 
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721  CONTINUE 
C 

DO  36S  J=  1.256.2 

\VRITE(24,23)  J.INT(129.J).J+l.INT(r29.J+l) 
23  F  ORNIAT  (2X.2(I3,3X,  1PE 1 1  .4.3X)) 

368  CONTINUE 
C 

C  FIND  THE  POINT  OF  MAXIMUM  INTENSITY 
C 

M\X=0.0 
DO  124  1=1.256 
DO  123  J= 1.256 
IF(INT(I.J).GT.M\X)  THEN 
IM\X=I 
JMAX=J 
\L\X=INT(I,J) 

END  IF 

123  CONTINUE 

124  CONTINUE 
C 

C 

C  SELECT  BOX  TO  GRAPH 
C 

IOFFSET =DL\X-48-l 
J  OFFSE  T=JN  L-VX-48- 1 
C 

IF(IOFFSET.LT.l)  THEN 
IOFFSET=0 
END  IF 

IF(IOFFSET+96.GT.256)  THEN 
IOFFSET=256-96-l 
ENDIF 
C 

IF(JOFFSET.LT.l)  THEN 
JOFFSET=0 
ENDIF 

IF(JOFFSET+96.GT.256)  THEN 
JOFFSET=256-96-i 
ENDIF 
C 

CM.-VX=0.0 
DO  542  1=1,96 
DO  541  J = 1 ,96 

M(I,J)=INT(IOFFSET+I.JOFFSET+J) 
IF(M(I,J).GT.CMVX)  THEN 
\1AXI=I 
MAXJ=J 
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C.\iAX=M(I.J) 
END  IF 
CONTINUE 
CO  NT  I  NIT 


OUTPUT  OF  INTENSITY 

\VRITE(26,73)  INIAX.JMAX.MAX 
FORMAT(2X,I3,3X.I3.3X,lPEll.4) 
DO  702  1=1.96 
DO  701  J=  1.96 
WRITE(26.73)  I.J.M(I.J) 
CONTINUE 
CONTI  NET 
STOP 
END 


APPENDIX  C 


d 


This  appendix  contains  the  code  for  the  program  ZERNIKF.  This  program 
calculates  the  far- field  intensity  pattern  of  the  impulse  response  for  the  two  lenses 
at  various  angles  of  incident  radiation.  The  aberration  coefficients  were  first  cal¬ 
culated  by  the  ray  trace  program  FALCON,  and  then  these  coefficients  were 
mapped  onto  a  unit  circle  which  is  the  exit  pupil  of  the  optical  system.  These 
values  were  then  entered  into  a  Fast  Fourier  Transform,  and  the  intensity  plots 
were  calculated  by  squaring  the  output  of  the  transform.  The  optic  axis  is  cen¬ 
tered  on  the  position  (128.5,128.5)  in  the  object  plane,  and  in  the  image  plane  the 
system  rounds  the  optic  axis  to  be  at  the  position  (129.129). 
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P R OG RAM:  \ ! Z E R XI KE . F 


C 
C’ 

C' 

c 

C  THIS  PROGR.VM  CALC  ULATES  THE  INTENSITY  OF  THE  TWO 

C  DIMENSION. VI.  IMAGE  OF  A  PLANE  WAVE  THROUGH  A 

C  LEX'S  WITH  THE  .VPPROPRLYTE  ZERX1KE  ABERRATIONS  AT  SOME 
ANGLE  ALPHA 


I  NT  EGER  I A 1 . 1 A  AN  I .  X2 .  N3. 1 J  OB.IWK(2000 ) .  N 
IXTEGE R  L  YR Ci.JAR G.CENTE R 
RE.AL  R\VKi20O0).LXT(256.256).X.Y.M(256.256).R  1 .112 
REAL  IXDEX.PI.RAD ARG2.FOC.YL.M-YX.il. LAMBDA 
COMPLEX  C'WK(256),A(256,256).B(256.256) 

OPEN(UNIT=23.$TATUS=’NEW\FILE=’APER.l)AT‘  > 
0PEN(UNIT=2LSTATUS=,NEW\FILE='DBGEET.I)AT) 
OPEN(UNIT=25,STATUS=ENEWEFILE='TEST.I)AT) 
OPEN(UNIT=26,STATUS— 'NEW' .FILE='ZER  NIKE.  DAT) 

RE.AD  PARAMETERS 

N=256 

LAMBDA =.6328 
PI=3. 1415927 
\VRITE(5,29) 

RE.AD(5,30)  Rl 
RE.AD(5.30)  R2 
READ(5,30)  .ALPELA 
READ(5.30)  INDEX 
READ(5,30)  A1 
RE.AD(5,30)  A2 
RE  AD  (5, 30)  A3 
READ(5,30)  A-l 
READ(5,30)  Ao 
READ(5,30)  A 6 
READ(o,30)  A 7 
RE  AD  (5,30)  A8 
RE.AD(5,31)  CENTER 

29  FORNL\T(2X,'TYTE  IN  YAJT'ES  FOR  RLR2.ALPHA.INDEX.A1-8.  AND 

*  CENTER(0  OR  1)  RETURN  AFTER  EACH  ENTRY,  REMEMBER  PERIODS!') 

30  FOR.\L\T(F  15.-1) 

31  FOR\LYT(l3) 

WRITE(5,30)  Rl 

C 

C  Rl=RADIUS  OF  CURVATURE  OF  LEFT  HAND  SURFACE 
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C  R2=RAD  I  US  OF  CTPYATl'RE  OF  RIGHT  ILVXD  SIP, FACE 

C  1XDEX=LXDEX  OF  REFRACTION  OF  THE  LENS  MEDIUM 
C  F OC.AL= FOCAL  LENGTH  OF  THE  THIN  LENS 

C  A1-S=THE  ZERNIKE  COEFFICIENTS  :  TILT.FOC'US.ASTIG\L\TISM.CONL\. 
C  AND  SPHERICAL 

C  CENTER= WHETHER  OR  NOT  YOU  WANT  TITE  INTENSITY  GRAPH  TB 
C  CENTERED  ON  THE  NLAXLMUM  OR  NOT 

C 

F  OCAL=  1  ./((INDEX- 1  .)*(  I /Rl- 1  /R2)) 

ALPHA=ALPILV*PI/130. 

C 

C  CREATE  .APERT LPE 
C 

DO  220  I=I.N 
X3=I 

DO  200  J=1,N 
YJ=J 

RAD=SQRT((XI-(l27.-F0.5))**2+(YJ-(127.-f0.5))**2)/17. 

X=(XI-(  127.4-0. 5))/ 17. 

Y=(YJ-(127.4-0.5))/l7. 

IF(RAD.LT.l.)  GO  TO  60 
IF(RAD.GT.l.)  GO  TO  100 

60  PHI  =■  A1*Y  +  A2*X  +  A3*(-1.+2.*Y**2+2.*X**2) 

*  +  A5*2.*X*Y  +  A4*(Y**2-X**2) 

*  +  A6*(-2.*Y+3.*Y**3+3.*Y*X**2) 

*  +  A7*(-2.*X+3.*X*Y**2+3.*X**3) 

*  +  A8*(1.-6.*Y**2-6.*X**2+6.*Y**4+ 

*  12.*X**2*Y**2+6.*X**4) 

ARG2=2.*PI*PHI  -  2.*PI*Y*SIN(JALPH-\)*10**4/L.AMBDA 
A(I,J)=CMPLX(COS(ARG2).SIN(ARG2)) 

GO  TO  120 

100  A(I,J)=(0. 0.0.0) 

120  CONTINUE 

200  CONTINUE 
220  CONTINUE 
C 

C  OUTPUT  APERTURE  MATRIX 
C 

DO  230  1=90.127 
DO  225  J=90.127 
WRITE(23.40)  L.I.A(L.J) 

10  FORMAT(2X.l3.3X.i3.2(3X.  El'll  1  t 

225  CONTINUE 
230  CONTINUE 
C 
C 

C  PERFORM  F  \-T  HO."  ■ 


V 

, 

1 

! 

i 

c 

' 

c 

1 

240 

1 

L\1=N 

IA2=N 

N1=N 

.  > 

N2=N 

$ 

N3=l 

% 

.f 

1 

1  C 

IJOB=l 

CALL  FFT3D  (A,LAl,IA2,Nl,N2,N3,IJOB,I\VK,R\VK,C\VK) 

>*■ 

1 

i  c 

OUTPUT  NEW  MATRIX  A 

!• 

*N 

c 

DO  242  J=1,N 

■ 

1=127 

WRITE(24,40)  I,J,A(I,J) 

V 

5 

J  40 

F0RMAT(2X,I3,3X,I3,2(3X,1PE11.4)) 

V 

♦a 

7  242 

CONTINUE 

V 

c 

\ 

c 

C 

EF(CENTER.EQ.l)  GO  TO  767 

i« 

■j 

IF(CENTER.EO.O)  THEN 

J 

5  c 

c 

SHIFT  ZERO  FREQUENCY  TO  ARRAY  CENTER 

c 

c 

J 

DO  260  1=1,256 

* 

■%. 

\ 

DO  255  J=  1 ,128 

n 

&1 

i 

B(I,J)=A(I,J+128) 

>* 

255 

B(I,J+128)=A(I,J) 

CONTINUE 

260 

CONTINUE 

DO  270  J=l,256 

■> 

5 

8 

DO  265  1=1,128 

[> 

i 

C 

c 

A(I,J)=B(I+128,J) 

A(I+128,J)=B(I,J) 

265 

CONTINUE 

<** 

270 

CONTINUE 

> 

r 

i  c 

*■  t 

l  c 

c 

c 

REDUCE  SIZE  OF  ARRAY 

c 

V 

c 

DO  250  1=1,96 

i 

j3 

DO  245  J=l,96 

c 

t 

u 

B(I,J)=A(I+80,J+80) 

1 

245 

CONTINUE 

250 

CONTINUE 

c 

a 

C 

88 

1* 

V 

-- 

l4 

i*1 

*'  ♦*.  **•:*•  >  >'  •  .  •>  **.  *».y.  t  <■  -  ■ 

c 

C  WRITE  CENTERED  A  MATRIX 
C 

DO  274  1=1,96 
DO  272  J=l,96 
WRITE(25,40)  I,J,B(I,J) 

272  CONTINUE 
274  CONTINUE 
C 
C 

C  CALCULATE  THE  MODULUS  AND  INTENSITY  OF  ARRAY 

C 

C 

DO  400  1=1,96 
DO  380  J=l,96 
M(I,J)=CABS(B(I,J)) 

INT(I,J)=M(I,J)**2 
380  CONTINUE 
400  CONTINUE 
C 
C 

C  FIND  THE  POINT  OF  MAXIMUM  INTENSITY 
C 

MAX=0.0 
DO  343  1=1,96 
DO  342  J=l,96 
1F(INT(I,J).GT.MAX)  THEN 
IMAX=I 
JMAX=J 
MAX=INT(I,J) 

END  IF 

342  CONTINUE 

343  CONTINUE 
C 

C 

C  OUTPUT  OF  INTENSITY 
C 

IMAX=LMAX+80 
JMAX=  JMAX+SO 
WRITE (26,134)  I\LAX,JMAX,MAX 
134  FORMAT(2X,I3,3X,I3,3X,lPEll.4) 

DO  404  1=1,96 
DO  402  J=l,96 
WRITE(26,41)  I,J,INT(I,J) 

4 1  FORMAT(2X.I3,3X,I3,3X,lPEl  1 .4) 

402  CONTINUE 


404  CONTINUE 

C 

C 

END  IF 
GO  TO  999 

Q  ************************************************************************ 

C 

C  TO  GRAPH  MAXIMUM  INTENSITY  AT  CENTER  OF  PLOT 
C 

C  SHIFT  ZERO  FREQUENCY  TO  ARRAY  CENTER 
C 

767  DO  654  1=1,256 
DO  653  J=l,128 
B(I,J)=A(I,J+128) 

B(I,J+128)=A(I,J) 

653  CONTINUE 

654  CONTINUE 
DO  652  J=l,256 

DO  651  1=1,128 
A(I,J)=B(I+128,J) 

A(I+128,J)=B(I,J) 

651  CONTINUE 

652  CONTINUE 
C 

C 

C  CALCULATE  THE  MODULUS  AND  INTENSITY  OF  THE  ARRAY 

C 

C 

DO  721  1=1,256 
DO  720  J= 1,256 
M(I,J)=CABS(A(I,J)) 

INT(I,J)=M(I,J)**2 

720  CONTINUE 

721  CONTINUE 
C 

C 

C  FIND  THE  POINT  OF  MAXIMUM  INTENSITY 
C 

NL\X=0.0 
DO  124  1=1,256 
DO  123  J=  1.256 
IF(INT(I,J).GT.MAX)  THEN 
I\LAX=I 
JMAX=J 
MAX=INT(I,J) 

ENDIF 

123  CONTINUE 
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124  CONTINUE 

C 

C 

C  SELECT  BON  TO  GRAPH 
C 

IOFFSET=IMAX-48-l 
J  OFF  SET= JMAX-48- 1 
C 

IF(IOFFSET.LT.l)  THEN 
IOFFSET=0 
ENDIF 

IF(IOFFSET+96.GT.256)  THEN 
IOFFSET=256-96-l 
ENDIF 
C 

IF(JOFFSET.LT.l)  THEN 
JOFFSET=0 
ENDIF 

IF(JOFFSET+96.GT.256)  THEN 
JOFFSET— 256-96-1 
ENDIF 
C 

CMAX=0.0 
DO  542  1=1,96 
DO  541  J=l,96 

M(I,  J)=INT(IOFFSET +1,  JOFFSET + J) 
IF(M(I,  J).GT.CMAX)  THEN 
MAXI = I 
MAXJ=J 
CMAX=M(I,J) 

ENDIF 

541  CONTINUE 

542  CONTINUE 
C 

C 

C  OUTPUT  OF  INTENSITY 
C 

WRITE(26,34)  IMAX,JMAX,MAX 
34  FORNL\T(2X,I3,3X,I3,3X,lPEll.4) 

DO  702  1=1,96 
DO  701  J=l,96 
WRITE(26,73)  I,J,M(I,J) 

73  F0RMAT(2X,I3,3X,I3,3X,1PE1 1.4) 

701  CONTINUE 

702  CONTINUE 
999  STOP 

END 
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APPENDIX  D 


This  appendix  contains  a  representative  code  for  multi-aperture  systems. 
The  program  given  is  MULTI.  This  program  computes  the  multi-aperture 
impulse  response  for  the  Walton  thin  lens  at  various  angles  a  for  various  numbers 
of  apertures  at  various  placements  in  the  field.  Each  lens  or  aperture  is  assumed 
to  be  exactly  like  any  other  aperture  in  the  system.  The  aperture  knowns  are  the 
same  as  in  the  program  ALPHA.  After  the  program  has  calculated  the  appropri¬ 
ate  distances  and  therefore  phase  lags  for  one  aperture,  it  then  translates  this 
aperture  to  various  determined  positions  in  the  exit  plane  of  the  optical  system. 
After  the  multi-aperture  system  is  set  up  in  the  object  plane,  the  field  is  then 
Fresnel  propagated  to  the  image  plane  where  the  intensity  plots  are  computed, 
and  given  in  the  text.  The  parameters,  as  well  as  other  dimensions  determined  by 
the  computer,  are  in  units  of  pixels.  The  optic  axis  is  centered  on  the  position 
(128.5,128.5)  in  the  object  plane,  and  in  the  image  plane  the  system  rounds  the 
optic  axis  to  be  at  the  position  (129,129). 
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PROGRAM:  MULTI.F 


THIS  PROGRAM  CALCULATES  THE  INTENSITY  OF  THE  TWO 
DIMENSIONAL  IMAGE  OF  A  PLANE  WAVE  THROUGH  A 
MULTI  APERTURE  SYSTEM  USING  THE  WALTON  THIN 
LENS  AT  SOME  ANGLE  ALPHA 


INTEGER  IAlXA2,Nl,N2,N3,IJOB,IWK(2000),N 

INTEGER  LARG,JARG, FI, F2,APE, CENTER 

REAL  RWK(2000),INT(256,256),X,Y,M(256,256),R1,R2,H,ALPHA 

REAL  INDEX.PI, RAD, BETA, ETA, XSIM,XSIP,D,ARG1,XRG2F, ROT 

REAL  LAMBDA WS TOP ,MAX,ARG3 

COMPLEX  CWK(256),A(256, 256)3(256,256) 

OPEN(UNIT=23,STATUS=’NEW’,FILE=’APER.DAT’) 

OPEN(UNIT=24,STATUS=’NEW’,FILE=’DBGFFTX)AT') 

OPEN(UNIT=25,STATUS=’NEW’,FILE=’MULTI.DAT’) 

OPEN(UNIT=26,STATUS=’NEW\FILE=,MULTI.DAT’) 

OPEN(UNIT=27,STATUS=’NEW’,FILE=’MILLSl.DAT’) 

READ  PARAMETERS 


N=256 
PI=3. 1415927 
\VRITE(5,29) 

READ(5,30)  Rl 
READ(5,30)  R2 
READ(5,30)  H 
READ(5,30)  ASTOP 
READ(5,30)  R.ADIUS 
READ(5,30)  APER 
READ(5,30)  ROT 
READ(5,30)  ALPHA. 

READ(5,30)  INDEX 
READ(5,31)  CENTER 

29  FORNLA.T(2X,'TA'PE  IN  V.ALUES  FOR  RlR2,H,  ASTOP  .RAD  IUS.APERROII 
*  INDEX, AND  CENTER,  RETURN  .AFTER  EACH  ENTRY,  REMEMBER  PERIODS!') 

30  F0RMAT(F15.4) 

31  F0RN1AT(I3) 

C 

C  R1=R.ADIUS  OF  CURVATLRE  OF  LEFT  HAND  SERF  ACE 

C  R2=RADIUS  OF  CURVATURE  OF  RIGHT  HAND  SURFACE 

C  H=HEIGHT  OF  LENS  FROM  OPTIC  .AXIS  AT  CENTER  OF  LENS 
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C  ASTOP=RADIUS  OF  THE  APERTURE  STOP  OF  THE  SYSTEM 
C  RADrUS=OUTER  RADIUS  OF  SYNTHESIZED  APERTURE 
C  (EXPRESSED  AS  MULTIPLE  OF  SUB-APERTURE  RADIUS) 

C  APER = NUMBER  OF  SUB- APERTURES  IN  THE  ARRAY 

C  ROT=ROTATION  OF  ARRAY  (DEGREES) 

C  ALPHA= ANGLE  WITH  RESPECT  TO  OPTIC  AXIS  OF  INCIDENT 

C  PLANE  WAVE 

C  INDEX=INDEX  OF  REFRACTION  OF  THE  LENS  MEDIUM 

C  LAMBDA*  W'AVELENGTH  OF  INTEREST 
C  F=FOCAL  LENGTH  OF  THE  THIN  LENS 
C 

ALPHA*  ALPHA*PI/ ISO. 

DELX=ASTOP/17. 

LAMBDA*. 6328 
ROT=ROT*PI/1SO. 

F=1./((INDEX-1.)*(1/R1-1/R2)) 

SDIAM=2.*17. 

APE=APER-1. 

ORAD*RADIUS*17. 

WRITE(5,30)  ALPHA 

CREATE  APERTURE 


FIRST  APERTURE 

IARG=128-17.+(ORAD-17.)*COS(ROT) 
JARG*128-17.+(ORAD-17.)*SIN(ROT) 

DO  107  1=1, N 
DO  105  J=1,N 
A(I,J)=CMPLX(0,0) 

105  CONTINUE 
107  CONTINUE 
C 

DO  220  I=1,SDIAM 
XI = I 

DO  200  J=1,SDIAM 
YJ=J 

RAD=DELX*SQRT((XI-(17.-H.5))**2.+(YJ-(17.-H.5))**2.) 
X=(X3-(17.+1.5))*DELX 
Y=(YJ-(17.+1.5))*DELX 
IF(RAD.LT.ASTOP)  GO  TO  60 
IF(RAD.GT.ASTOP)  GO  TO  100 
60  BETA=l./(COS(ALPR\)**2.) 

ETA=X*  *2.+Y*  *2.-H*  *2. 
XSIM=(SQRT(R2**2.-H**2.))-Y*TAN(ALPR\) 
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XSIP=(SQRT(Rl**2.-H**2.))+Y*TAN(  ALPHA) 
D=COS(ALPHA)*((SQRT(XSIM**2.-ETA*BETA))+ 

i 

(S  QR  T(XSIP  *  *  2  .-E  T  A  *  BE  T  A))-XSIM-XS  IP ) 
ARG1=PI*(X**2.+Y**2.)/F 

ARG2=2.*PI*(INDEX-1.)*D+ARG1-2.*PI*X*SIN(ALPHA) 
ARG3= ARG2*  10**4 /LAMBDA 

IF(I.EQ.127)  THEN 

\VRITE(27,98)  J,BETA,ETA,XSIM,XSIP,D,ARGl,ARG2 

98 

1 

F0RMAT(2X,I3,7(3X,1PE11.4)) 

ENDIF 

A(I+IARG,J+JARG)=CMPLX(COS(ARG3),SIN(ARG3)) 

GO  TO  120 

100 

A(I+LARG,J+JARG)=(0.0,0.0) 

120 

CONTINUE 

m  200 

CONTINUE 

ft  220 

CONTINUE 

I  C 

CREATE  REMAINING  SUB-APERTURES 

3  c 

« 

DO  227  1=1, APE 

:  is 

U=(ORAD-17.)*COS(I*2.*PI/APER+ROT) 

V=(ORAD-17.)*SIN(I*2.*PI/APER+ROT) 

XP=127.-17.+U 

YP=127.-17.+V 

DO  225  J=1,SDIAM 

V. 

DO  223  K=1,SDIAM 

A(XP+J,YP+K)=A(XP+J,YP+K)+A(J+LVRG,K-f- JARG) 

223 

CONTINUE 

225 

CONTINUE 

H  -»  227 

i  S  c 

CONTINUE 

I  C 

OUTPUT  APERTURE  MATRIX 

1  as 

DO  230  1=110,127 

DO  229  J=1 10,127 

i  1 

\VRITE(23,40)  I,J,A(I,J) 

4  40 

F0RMAT(2X,I3,3X,I3,2(3X,1PE11.4)) 

229 

CONTINUE 

230 

P 

CONTINUE 

I  c 

c  c 

PERFORM  FAST  FOURIER  TRANSFORM  USING  IMSL 

c 

c 
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240  IA1=N 
IA2=N 
Nl=N 
N2=N 
N3=l 
IJOB=l 

CALL  FFT3D  (A,IAl,IA2,Nl,N2,N3,IJOB,IWK,RWK,CWK) 
C 

C  OUTPUT  NEW  MATRIX  A 
C 

DO  242  J=1,N 
1=127 

WRITE(24,40)  I,J,A(I,J) 

40  FORMAT(2X,I3,3X,I3,2(3X,lPEll.4)) 

242  CONTINUE 

C 

C 

IF(CENTER.EQ.l)  GO  TO  767 
IF(CENTER.EQ.O)  THEN 
C 
C 

C  SHIFT  ZERO  FREQUENCY  TO  ARRAY  CENTER 

C 

C 

DO  260  1=1,256 
DO  255  J=l,128 
B(I,J)=A(I,J+128) 

B(I,J+128)=A(I,J) 

255  CONTINUE 
260  CONTINUE 
DO  270  J=  1,256 
DO  265  1=1,128 
A(I,J)=B(I+128,J) 

A(I+128,J)=B(I,J) 

265  CONTINUE 
270  CONTINUE 
C 
C 

C  REDUCE  SIZE  OF  ARRAY 

C 

C 

DO  250  1=1,96 
DO  245  J=l,96 
B(I,J)=A(I+80,J+80) 

245  CONTINUE 
250  CONTINUE 
C 
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c 

C  WRITE  CENTERED  A  MATRIX 
C 

DO  274  1=1,96 
DO  272  J=l,96 
WRITE(25,40)  I,J,B(I,J) 

272  CONTINUE 
274  CONTINUE 
C 
C 

C  CALCULATE  THE  MODULUS  AND  INTENSITY  OF  ARRAY 

C 

C 

DO  400  1=1,96 
DO  380  J=l,96 
M(I,J)=CABS(B(I,J)) 

INT(I,J)=M(I,J)**2 
380  CONTINUE 
400  CONTINUE 
C 
C 

C  FIND  THE  POINT  OF  MAXIMUM  INTENSITY 
C 

MAX=0.0 
DO  343  1=1,96 
DO  342  J=l,96 
EF(INT(I,J).GT.MAX)  THEN 
IMAX=I 
JMAX=J 
MAX=INT(I,J) 

END  IF 

342  CONTINUE 

343  CONTINUE 
C 

C 

C  OUTPUT  OF  INTENSITY 
C 

IMAX=MAX+80 
JMAX=JMAX+80 
WRITE(26,134)  IM\X,JMAX,MAX 
134  F0RM\T(2X,I3,3X,I3,3X,1PE11.4) 

DO  404  1=1,96 
DO  402  J=l,96 
WRITE(26,41)  I,J,INT(I,J) 

4 1  F0RM\T(2X,I3,3X,I3,3X,  1  PEI  1 .4) 

402  CONTINUE 
404  CONTINUE 
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i 


I 


i 


i 

E 


c 

END  IF 
GO  TO  999 
C 

Q  ***************************************************************** 

c 

C  TO  GRAPH  MAXIMUM  INTENSITY  AT  CENTER  OF  PLOT 
C 

C  SHIFT  ZERO  FREQUENCY  TO  .ARRAY  CENTER 

C 

C 

767  DO  654  1=1,256 
DO  653  J=  1,128 
B(I,J)=A(I,J+128) 

B(I,J+128)=A(I,J) 

653  CONTINUE 

654  CONTINUE 
DO  652  J=  1.256 

DO  651  1=1,128 
A(I,J)=B(I+128,J) 

A(I+128,J)=B(I,J) 

651  CONTINUE 

652  CONTINUE 
C 

C 

C  CALCULATE  THE  MODULUS  AND  INTENSITY  OF  THE  .ARRAY 

C 

C 

DO  721  1=1,256 
DO  720  J=  1,256 
M(I,J)=CABS(A(I,J)) 

INT(I,J)=M(I,J)**2 

720  CONTINUE 

721  CONTINUE 
C 

C 

C  FIND  THE  POINT  OF  MAXIMUM  INTENSITY 
C 

MAX=0.0 
DO  124  1=1,256 
DO  123  J=  1,256 
IF(INT(I,J).GT.MAX)  THEN 
I  MAX = I 
J\LAX=J 
MAX=INT(I,J) 

ENDIF 

123  CONTINU'D 


c 
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124  CONTINUE 

C 

C 

C  SELECT  BOX  TO  GRAPH 
C 

IOFFSET =I\LXX-48-l 
J  OFF  SET= JMAX-48- 1 
C 

IF(IOFFSET.LT.l)  THEN 
IOFFSET=0 
ENDIF 

IF(IOFFSET+96.GT.256)  THEN 
IOFFSET =256-96-1 
ENDIF 
C 

IF(JOFFSET.LT.l)  THEN 
JOFFSET=0 
ENDIF 

IF(JOFFSET+96.GT.256)  THEN 
JOFFSET =256-96- 1 
ENDIF 
C 

DO  542  1=1,96 
DO  541  J=l,96 

M(I,J)=INT(IOFFSET+I,JOFFSET+J) 

541  CONTINUE 

542  CONTINUE 
C 

C 

C  OUTPUT  OF  INTENSITY 
C 

WRITE(26,73)  IMAX,JMAX,MAX 
73  F0RMAT(2X,I3,3X,I3,3X,1PE11.4) 

DO  702  1=1,96 
DO  701  J=l,96 
WRITE(26,73)  I,J,M(I,J) 

701  CONTINUE 

702  CONTINUE 
999  STOP 

END 
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